MIT-CTP-3807 



The rigidity of crystalline color superconducting quark matter 

Massimo Mannarelli,^'EI Krishna Rajagopal,^'[l] and Rishi Sharma^'0 

' Center for Theoretical Physics, Massachusetts Institute of Technology, Cambridge, MA 02139 

(Dated: August 7, 2007) 

We calculate the shear modulus of crystalline color superconducting quark matter, showing that 
this phase of dense, but not asymptotically dense, three-flavor quark matter responds to shear stress 
like a very rigid solid. To evaluate the shear modulus, we derive the low energy effective Lagrangian 
that describes the phonons that originate from the spontaneous breaking of translation invariance 
by the spatial modulation of the gap parameter A. These massless bosons describe space- and 
time-dependent fluctuations of the crystal structure and are analogous to the phonons in ordinary 
crystals. The coefficients of the spatial derivative terms of the phonon effective Lagrangian are 
related to the elastic moduli of the crystal; the coefflcients that encode the linear response of the 
crystal to a shearing stress define the shear modulus. We analyze the two particular crystal structures 
which are energetically favored over a wide range of densities, in each case evaluating the phonon 
effective action and the shear modulus up to order in a Ginzburg-Landau expansion, finding 
shear moduli which are 20 to 1000 times larger than those of neutron star crusts. The crystalline 
color superconducting phase has long been known to be a superfluid — by picking a phase its order 
parameter breaks the quark-number U{1)b symmetry spontaneously. Our results demonstrate that 
this superfluid phase of matter is at the same time a rigid solid. We close with a rough estimate 
of the pinning force on the rotational vortices which would be formed embedded within this rigid 
superfluid upon rotation. Our results raise the possibility that (some) pulsar glitches could originate 
within a quark matter core deep within a neutron star. 
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I. INTRODUCTION 

Quantum chromodynamics (QCD) predicts that 
if one can squeeze nuclear matter to the point that 
the number of quarks per unit volume is greater than 
that within a nucleon, the matter that results cannot 
be described as made up of nucleons: the individual 
nucleons are crushed into quark matter It re- 
mains an open question whether "liberated" quark 
matter exists at the core of neutron stars, which 
have central densities at most ~ 10 times that of 
ordinary nuclear matter. One route to answering 
this question would be an ab initio calculation of 
the range of densities at which the expected first or- 
der transition from nuclear to quark matter occurs 
in QCD. No currently available calculational method 
promises such a result any time soon, which is not 
surprising since locating such a first order transition 
requires quantitative control of the pressure of both 
phases, making this a much more difficult question 
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than understanding the properties of either phase 
on its own. The more promising route to answering 
whether neutron stars feature quark matter cores re- 
quires first understanding the qualitative and semi- 
quantitative properties of quark matter and then dis- 
cerning how its presence deep within a neutron star 
would affect astronomically observed or observable 
phenomena. 

Except during the first few seconds after their 
birth, the temperature of neutron stars is expected 
to be of the order of tens of keV at most. At such 
low temperatures, if deconfined quark matter ex- 
ists in the core of neutron stars it is likely to be 
in one of the possible color superconducting phases, 
whose critical temperatures are generically of order 
tens of MeV 0. Furthermore, since compact star 
temperatures are so much smaller than these crit- 
ical temperatures, for many purposes (including in 
particular for a calculation of its rigidity) the quark 
matter that may be found within neutron stars is 
well-approximated as having zero temperature, as 
we shall assume throughout. 

The essence of color superconductivity is quark 
pairing, driven by the BCS mechanism which oper- 
ates whenever there is an attractive interaction be- 
tween fermions at a Fermi surface [3J. The QCD 
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quark-quark interaction is strong and attractive be- 
tween quarks that are antisymmetric in color, so we 
expect cold dense quark matter to generically exhibit 
color superconductivity. It is by now well established 
that at asymptotic densities, where the up, down 
and strange quarks can be treated on an equal foot- 
ing and the disruptive effects of the strange quark 
mass can be neglected, quark matter is in the color 
flavor locked (CFL) phase 0,11]. The CFL conden- 
sate is antisymmetric in color and flavor indices and 
therefore involves pairing between quarks of different 
colors and flavors. 

To describe quark matter as may exist in the cores 
of compact stars, we need to consider quark number 
chemical potentials /i of order 500 MeV at most (cor- 
responding to fiB = 3/^ of order 1500 MeV). The rel- 
evant range of /i is low enough that the strange quark 
mass, Ms, lying somewhere between its current mass 
of order 100 MeV and its vacuum constituent mass 
of order 500 MeV, cannot be neglected. Further- 
more, bulk matter, if present in neutron stars, must 
be in weak equilibrium and must be electrically and 
color neutral. All these factors work to separate the 
Fermi momenta of the quarks by ~ / fj,, and thus 
disfavor the cross-species BCS pairing which charac- 
terizes the CFL phase. If we imagine beginning at 
asymptotically high densities and reducing the den- 
sity, and suppose that CFL pairing is disrupted by 
the heaviness of the strange quark before color su- 
perconducting quark matter is superseded by bary- 
onic matter, the CFL phase must be replaced by 
some other superconducting phase of quark matter 
which has less, and less symmetric, pairing. 

Considering only homogeneous condensates, the 
next phase down in density is the gapless CFL 
(gCFL) phase @, i, 0, i, i, M MMM- In 
this phase, quarks of all three colors and all three 
flavors still form ordinary Cooper pairs, but there 
are regions of momentum space in which certain 
quarks do not succeed in pairing, and these regions 
are bounded by momenta at which certain fermionic 
quasiparticles are gapless. This variation on BCS 
pairing — in which the same species of fermions that 
pair feature gapless quasiparticles — was first pro- 
posed for two flavor quark matter [T3 | and in an 
atomic physics context [l^. In all these contexts, 
however, the gapless paired state turns out in general 
to suffer from a "magnetic instability" : it can lower 
its energ y b y the formation of counter-propagating 
currents [TgI . [iTj . 

In Refs. [H, [l^, [2^ the possibility that the next 



superconducting phase down in density has a crys- 
talline structure was explored. (For other possible 
patterns of three-flavor and three-color pairing, see 
Refs. [21I, [22|. For possibilities in which only two 
colors of quarks pair, see Refs. [2^. For possibil- 
ities in which only quarks of the same flavor pair, 
see Refs. [23|. See Ref. [2^ for an exhaustive study 
of less symmetric phases with BCS pairing only, the 
simplest example being the 2SC phase [2y, [23|, in- 
dicating that none are favorable.) Crystalline color 
superconductivity — the QCD analogue of a form of 
non-BCS pairing first considered by Larkin, Ovchin- 
nikov, Fulde and Ferrell (LOFF) [H] — is an at- 
tractive candidate phase in the intermediate den- 
sity regime [1^, HO, [U, [H, |H, H^l, because it al- 
lows quarks living on split Fermi surfaces to pair 
with each other. It does so by allowing Cooper pairs 
with nonzero total momentum 2q'^, for some set of 
q'^'s whose magnitude is fixed and is of order the 
splitting between Fermi surfaces, and whose direc- 
tions must be determined. In position space, this 
corresponds to condensates that vary in space like 
J2a cxp(2iq'* • x), meaning that the q^'s are the re- 
ciprocal vectors which define the crystal structure 
of the condensate which is modulated periodically 
in space and therefore spontaneously breaks space 
translation invariance. These phases seem to be free 
from magnetic instability |35j . which is consistent 
with the result that many of them have free ener- 
gies that are lower than that of the (unstable) gCFL 
phase for wide ranges of parameter values [20] • 

Two crystal structures for three-flavor crystalline 
color superconducting quark matter have been iden- 
tified as having particularly low free energy [20| . In 
the first, called the 2Cube45z phase, the {ud) and 
(lis) condensates each form face-centered cubes in 
position space. (In momentum space, each has eight 
q° vectors pointing toward the corners of a cube, 
which is the Bravais lattice of a face-centered cube.) 
The {us) cube is rotated relative to the (ud) cube by 
45 degrees. In the second favorable crystal structure, 
called the CubeX phase, the (ud) and (us) conden- 
sates each describe pairing with four different wave 
vectors q°, with the eight q^'s together forming the 
Bravais lattice of a single face-centered cube. The 
lattice spacing for these crystal structures is many 
tens of fm, as we shall discuss further in Section II, 
where we specify these crystal structures precisely. 

These phases were identified upon making two 
crucial approximations [20| , which we shall also em- 
ploy. First, the QCD interaction between quarks 
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was replaced by a point-like Nambu Jona-Lasinio 
(NJL) interaction with the quantum numbers of sin- 
gle gluon exchange. This is a conservative approach 
to the study of crystalline color superconductivity, 
since a point-like interaction tends to disfavor crys- 
talline color superconductivity relative to a forward- 
scattering dominated interaction like single gluon 
exchange [S^. Second, the calculation employed a 
Ginzburg-Landau expansion in powers of the gap pa- 
rameter A even though the values of A found in the 
calculation are large enough that the control param- 
eter for this expansion, A/(Af^/8/i), is about 1/2, 
meaning that the expansion is pushed to or beyond 
the boundary of its regime of quantitative validity. 
Upon making these approximations, one or other of 
the two most favorable crystalline color supercon- 
ducting phases is favored (i.e. has lower free energy 
than unpaired quark matter or any spatially uniform 
paired phase including the CFL and gCFL phases) 
for 

2.9 Ao < — < 10.4 Ao . (1) 

Here, Aq is the gap parameter in the CFL phase with 
Ms = 0, which is expected to lie somewhere between 
10 and 100 MeV 0. The estimate (P) is valid upon 
assuming A < M^/8fj,, Aq < fi and < ^j,"^, 
which here means 10.4 Aq < ^. In reality, both Mg 
and Ao are functions of /i, and the regime P]) trans- 
lates into a window in fi. See Ref. [36| for an anal- 
ysis in which both Ms{fJ,) and Ao(^) are calculated 
self-consistently in an NJL model, confirming that 
the window in fj. corresponding to ([T]) is broad in- 
deed. Even though it may not be quantitatively re- 
liable, the breadth of the regime H]), corresponding 
to the remarkable robustness of the two most fa- 
vorable crystalline phases which can be understood 
on qualitative grounds [20l |. makes it plausible that 
quark matter at densities accessible in neutron stars, 
say with ~ 350 — 500 MeV, will be in a crystalline 
phase. (Unless Aq lies at the upper end of its allowed 
range, in which case quark matter at accessible den- 
sities will be in the CFL phase.) 

In this paper, we study the elastic properties of 
the CubeX and the 2Cubc45z crystalline phases of 
three-flavor quark matter. We evaluate the shear 
moduli of these two structures by computing, in the 
Ginzburg-Landau approximation, the effective La- 
grangian for the phonon modes which emerge due 
to the spontaneous breaking of translation invari- 
ancc by the crystalline condensates. (See Ref. [3l| 



for an analysis of phonons in 2-flavor crystalline 
color superconducting phases.) The shear modulus 
is related to the coefficients of the spatial deriva- 
tive terms that appear in the phonon effective La- 
grangian. We find that when these two crystalline 
phases are subject to shear stresses, they behave like 
rigid sohds with shear moduli between 20 and 1000 
times larger than those characteristic of conventional 
neutron star crusts. 

We evaluate the shear modulus to order A^ in 
the Ginzburg-Landau expansion; higher order terms 
are suppressed only by powers of [A/(Afg /8^)]^, 
which is of order 1/4 in the most favorable crys- 
talline phases, large enough that significant higher 
order corrections to the shear modulus can be ex- 
pected. However, as we shall discuss below all that 
is required from the calculation of the shear modulus 
in making the case that pulsar glitches could origi- 
nate in a crystalline quark matter neutron star core 
is confidence that the shear modulus is comparable 
to or greater than that of the conventional neutron 
star crust. Since we find shear moduli that meet 
this criterion by a factor of 20 to 1000, there is cur- 
rently little motivation for evaluating higher order 
corrections. 

The possibility that quark matter could occur in a 
solid phase has been raised previously by Xu [iO] . He 
and his collaborators have explored some astrophys- 
ical consequences of a speculation that the quarks 
themselves could somehow be arranged in a crys- 
talline lattice. The crystalline color superconducting 
phase is very different in character: being a super- 
fluid, the quarks are certainly not frozen in place. 
Instead, what forms a crystalline pattern is the mag- 
nitude of the pairing condensate. Although it was 
clear prior to the present work that in this phase 
translational invariance is broken just as in a crystal, 
given that this phase is at the same time a super- 
fluid it was not clear (at least to us) whether it was 
rigid. Here, we demonstrate by explicit calculation 
that this phase, which as wc have discussed is plau- 
sibly the only form of quark matter that arises at 
densities reached within neutron star cores, is rigid 
indeed. Its shear modulus is parametrically of order 
A^^^, which could have been guessed on dimensional 
grounds. The shear modulus is in no way suppressed 
relative to this natural scale, even though the crys- 
talline color superconducting phase is superfluid. 

Our paper is organized as follows. In Section |TT] 
we describe the crystalline condensates which act 
as a background upon which the phonons act as 
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small perturbations, fully specifying the CubeX and 
2Cube45z structures. In Section III, after fully spec- 
ifying the NJL model that we shall use we intro- 
duce small displacements of a general diquark con- 
densate which breaks translational symmetries. We 
write a general expression for the effective action 
describing these displacement fields, by integrating 
out the fermions in the system. The final result for 
the phonon effective action is given in Eq. (|79|) . We 
relate the coefficients of the terms in the effective ac- 
tion involving the spatial derivatives of the displace- 
ment fields to the shear modulus in Subsection IIV Al 
and then evaluate these coefficients for the CubeX 
and 2Cube45z crystals in Subsections IIV CI and IIV Dl 
respectively. We end with a discussion of our results 
and their consequences in Section V. 

We have moved two relevant consistency checks to 
the Appendices to maintain continuity. In Appendix 
|X]we show explicitly that the displacement fields are 
massless to all orders in the gap parameter, as they 
must be by Goldstone's theorem. In Appendix |B] 
we evaluate the effective action for a simple "crys- 
talline" structure involving just two flavors of quarks 
and pairing with only a single wave vector q. In this 
case, the calculation can be done without making 
an expansion in the gap parameter, A. We flnd that 
the results in the limit of small A are consistent with 
the Ginzburg-Landau calculation. 

What are the observable signatures of the pres- 
ence of quark matter in a crystalline color super- 
conducting phase within the core of a neutron star? 
Prior to our work, the first steps toward calculat- 
ing the cooling rate for such neutron stars were 
taken in Ref. [4l|. Because the crystalline phases 
leave some quarks at their respective Fermi sur- 
faces unpaired, it seems likely that their neutrino 
emissivity and heat capacity will be only quantita- 
tively smaller than those of unpaired quark mat- 
ter [42, not parametrically suppressed. This sug- 
gests that neutron stars with crystalline quark mat- 
ter cores will cool by the direct URCA reactions, i.e. 
more rapidly than in standard cooling scenarios [43| . 
However, many other possible core compositions can 
open phase space for direct URCA reactions, mak- 
ing it unlikely that this will lead to a distinctive phe- 
nomenology [4^ . In contrast, the enormous rigidity 
that we have identified makes the crystalline phases 
of quark matter unique among all forms of matter 
proposed as candidates for the composition of neu- 
tron star cores. This makes it particularly impor- 
tant to investigate its consequences. We close this 



Introduction and close Section V of our paper with 
speculations about the consequences of a rigid quark 
matter core for puls ar g litch phenomenology, follow- 
ing the lead of Ref. [30| ■ 

A spinning neutron star observed as a pulsar 
gradually spins down as it loses rotational energy 
to electromagnetic radiation. But, every once in 
a while the angular velocity at the crust of the 
star is observed to increase suddenly in a dra- 
matic event called a glitch. The standard explana- 
tion m, m, \M, m, m m m es, m, y, y, m 113 

for their occurrence requires the presence of a super- 
fluid in some region of the star which also features a 
rigid array of spatial inhomogeneities which can pin 
the rotational vortices in the rotating superfluid. In 
the standard explanation of pulsar glitches, these 
conditions are met in the inner crust of a neutron 
star which features a crystalline array of positively 
charged nuclei bathed in a neutron superfluid (and 
a neutralizing fluid of electrons). The angular mo- 
mentum of the rotating superfluid is proportional to 
the density of vortices, meaning that as the star as 
a whole slows the vortices "want" to move apart. 
As they are pinned to a rigid structure, they can- 
not. Hence, this superfluid component of the star 
is spinning faster than the rest of the star. After 
sufficient time, when the "tension" in the vortices in 
some region reaches a critical point, there is a sudden 
"avalanche" in which vortices unpin and rearrange, 
reducing the angular momentum of the superfluid. 
The rest of the star, including in particular the sur- 
face whose angular velocity is observed, speeds up. 
We see that this mechanism requires superfluidity 
coexisting with a structure that is rigid enough that 
it does not easily deform when vortices pinned to it 
seek to move, and an adequate pinning force which 
pins vortices to the rigid structure in the flrst place. 

The crystalline color superconducting phases are 
superfluids: their condensates all spontaneously 
break the U{1)b symmetry corresponding to quark 
number. We shall always write the condensates as 
real. This choice of overall phase breaks U{1)b, and 
gradients in this phase correspond to supercurrents. 
And yet, as we show at length in this paper, crys- 
talline color superconductors are rigid solids with 
large shear moduli. Supcrsolids are another exam- 
ple of rigid superfluids [37|, i3^, i3S] , and on this ba- 
sis alone one might consider referring to crystalline 
color superconducting phases as supersolid. How- 
ever, supcrsolids are rigid by virtue of the presence 
of a conventional lattice of atoms whereas in a crys- 
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talline color superconductor individual quarks are 
not fixed in place in a rigid structure — it is the 
crystalline pattern of the diquark condensate that 
is rigid, exhibiting a nonzero shear modulus. This 
makes crystalline color superconductors sufficiently 
different from supersolids that we shall refrain from 
borrowing the term. 

Although our focus in this paper will be on the 
rigidity of crystalline color superconductors, it is 
worth pausing to consider how a condensate that 
vanishes at nodal planes can nevertheless carry a 
supcrcurrcnt. In an ordinary supcrfluid with spa- 
tially homogeneous condensate (j), the ground state 
(with no current) can be taken to have (jj real and 
the excited state with (jj = pexjpi9{r) carries a cur- 
rent proportional to p^S/O. If we now consider the 
simple example of a crystalline color superconductor 
with a diquark condensate of the form 

A(r) = Acos(2gz) , (2) 

this can carry current in the x- and y-directions in 
the usual fashion but how can a current in the z- 
direction flow through the nodal planes at which 
|A(r)| vanishes? The answer is that in adddition 
to making the condensate complex, with a spatially 
varying phase, the presence of a current must de- 
form the magnitude of the condensate slightly so as 
to make it nonvanishing at the location of the (for- 
mer) nodal planes. A deformation that has all the 
appropriate properties is 

A(r) = A(cos(2gz) + ie sin(2gz)) . (3) 

Note that if we write this condensate as A(r) = 
p{z) eyiY>i9{z) it satisfies p^d6/dz = 2e(7A^, mean- 
ing that the current is z-indcpcndent as required by 
current conservation.) The deformation of the con- 
densate is small for small e as is the current, which 
is proportional to e. At the location of the former 
nodal planes, the magnitude of the condensate is 
also proportional to e. We shall leave the determina- 
tion of whether ([3]) is the lowest energy deformation 
of the condensate with current 2eq/\'^ to future 
work; the correct current carrying state must be sim- 
ilar, even if it docs not turn out to be given precisely 
by ([3]) . We shall also defer the extension of this anal- 
ysis to the realistic, and more complicated, CubeX 
and 2Cube45z crystal structures. And, finally, we 
shall leave the analysis of currents flowing in circles 
in crystalline condensates — i.e. the construction of 
vortices — to future work. 



It is reasonable to expect that the superfluid vor- 
tices that will result when crystalline color super- 
conducting phases are rotated will have lower free 
energy if they are centered along intersections of the 
nodal planes of the underlying crystal structure, i.e. 
along lines along which the condensate already van- 
ishes even in the absence of a rotational vortex. By 
virtue of being simultaneously superfluids and rigid 
solids, then, the crystalline phases of quark matter 
provide all the necessary conditions to be the locus 
in which (some) pulsar glitches originate. The shear 
modulus, which we calculate, is larger that that of a 
conventional neutron star crust by a factor of 20 to 
1000, meaning that the rigidity of crystalline quark 
matter within a neutron star core is more than suf- 
ficient for glitches to originate therein. In Section V 
we also provide a crude estimate of the pinning force 
on vortices within crystalline color superconduct- 
ing quark matter and find that it is comparable to 
the corresponding value for neutron vortices within 
a neutron star crust. Together, our calculation of 
the shear modulus and estimate of the pinning force 
make the core glitch scenario worthy of quantita- 
tive investigation. The central questions that re- 
main to be addressed are the explicit construction 
of vortices in the crystalline phase and the calcula- 
tion of their pinning force, as well as the timescale 
over which sudden changes in the angular momen- 
tum of the core are communicated to the (observed) 
surface, presumably either via the common electron 
fluid or via magnetic stresses. 



II. THREE FLAVOR CRYSTALLINE 
COLOR SUPERCONDUCTIVITY 

A. Neutral unpaired quark matter 

We shall consider quark matter containing mass- 
less u and d quarks and s quarks with an effective 
mass Ms- The Lagrangian density describing this 
system in the absence of interactions is given by 

£o = M^f + pf 7o) , (4) 

where i,j — 1,2,3 are flavor indices and a,/3 = 
1,2,3 are color indices and we have suppressed the 
Dirac indices. The mass matrix is given by Mff^ = 

diag(0, 0, Ms)ij, whereas d^f = S'^f^Si^d and the 
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quark chemical potential matrix is given by 



Q/3 



(5) 

with Qij = diag(2/3, — 1/3, — l/3)y the quark 
electric-charge matrix and T3 and Ts the Gell-Mann 
matrices in color space, /i is the quark number chem- 
ical potential and if quark matter exists in neutron 
stars, we expect /i to be in the range 350 — 500 MeV. 

In QCD, /ie, /X3 and /Ltg are the zeroth components 
of electromagnetic and color gauge fields, and the 
gauge field dynamics ensure that they take on values 
such that the matter is neutral [Sq , satisfying 



0, 



(6) 



with Q the free energy density of the system. In the 
present paper we shall employ an NJL model with 
four-fermion interactions and no gauge fields. We in- 
troduce ^J■3 and Ms by hand, and choose them to 
satisfy the neutrality constraints ©. The assump- 
tion of weak equilibrium is built into the calculation 
via the fact that for a given color a the chemical 
potential of the d quark is equal the chemical poten- 
tial of the s quark which is equal to the sum of the 
chemical potential of the u quark and Me- 

Let us consider now the effect of a nonvanishing 
strange quark mass Ms- Suppose we start by set- 
ting the gauge chemical potentials to zero. Then, in 
weak equilibrium a nonzero Mg implies that there 
are fewer s quarks in the system than u and d quarks, 
and hence the system is positively charged. To re- 
store electrical neutrality, a positive Me is required 
which tends to reduce the number of up quarks rel- 
ative to the number of down and strange quarks. In 
the absence of pairing, color neutrality is obtained 
with M3 = Ms = 0. 

To lowest order in Af^ / fi^ , the effect of a nonzero 
strange quark mass can be taken into account by 
treating the strange quark as massless, but with a 
chemical potential that is lowered by M'^/{2^) from 
M + (Me/3). Indeed we can expand the Fermi mo- 
mentum of the strange quark as 



o 



M^ 



(7) 



= t;7, yielding 



/f 



Pf 



Pf 



Pf 



M H ^ = H 



M 

M - 

Ml 
4m 



Ml 
6m 



M2 

12m 4m 



(8) 



Now we need no longer to be careful about the dis- 
tinction between pp^s and m's, as we can simply 
think of the three flavors of quarks as if they have 
chemical potentials 



Md = Mu 

Mu = Pf 



2<5m2, 



with 



8m 



Sfi, 



(9) 



(10) 



and to this order electric neutrality requires that 



and we can write the chemical potential matrix as 
M^/ = ^"'^ ® diag (m« , Md , ^J■s ),j ■ ( 1 1 ) 

The factor 2 in front of the S^^s in Eq. ([9]) is taken 
to be consistent with the notation used in the anal- 
ysis of crystalline superconductivity in a two flavor 
model [29|, in which the two Fermi momenta were 
taken to be m i "^M meaning that they were sepa- 
rated by 2 Sfj,. In the three flavor case we have to 
contend with three combinations of pairs of Fermi 
surfaces and we define 2 Sfij = |Mi ^ Mfc| where j and 
k are flavor indices different from each other and 
from I. The subscripts on the Sfi^s in Eq. Q are 
consistent with this definition. The most important 
driver of qualitative changes in the physics of quark 
matter as a function of decreasing m is the increase 
in Sfj. = M'^/{8fi), the parameter which governs the 
splitting between the three Fermi momenta. 

Finally, we note that the equality of dfj,2 and Sfi^, 
is only valid to leading order in M^/m^- At the next 
order. Me = /(4m) ~ ^'^s^/(48m^) and hence Sfi^ = 
Me/2 while 6^2 = (^Ms + Mf/{16n^). The fact that 
(5M3 and 6fi2 arc close to, but not exactly, equal will 
be used in Appendix 1X1 in the explicit demonstration 
that the phonons are massless. The consequences of 
the fact that the splitting between the u and s Fermi 
surfaces is slightly larger than the splitting between 
the u and d Fermi surfaces were explored in Ref. [s^ . 
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B. The crystalline condensate 

The description in Section Hi Al is valid in the ab- 
sence of pairing between quarks. The interaction be- 
tween quarks is attractive between two quarks that 
are antisymmetric in color. This induces the for- 
mation of Cooper pairs 0, and the ground state 
develops a nonzero diquark condensate which is pre- 
dominantly antisymmetric in color indices. It is 
free-energetically favorable to have condensates that 
are antisymmetric in Dirac indices 0, 0, S HB] and 
consequently antisymmetric in flavor also. This im- 
plies that only quarks of different flavor form Cooper 
pairs. 

The fact that quarks that like to pair have Fermi 
surfaces that are split in the absence of pairing mo- 
tivates the possibility that quark matter may ex- 
ist as a crystalline superconductor in at least some 
part of the parameter space ,28, _29,, 2^ ■ This possi- 
bility was explored for three-flavor quark matter in 
Ref. [lOl (and for simpler condensates in [l^, [l9| ) 
where condensates of the form 

3 

1=1 q?e{qj 

were analyzed. As promised, the expression in 
Eq. ()12p is antisymmetric in color, flavor and Dirac 
indices. For a given / (which runs from 1 to 3), 
A/ represents the strength of the pairing between 
quarks whose color is not / and whose flavor is not 
/. The periodic modulation in space of the /'th con- 



densate is defined by a set of momentum vectors 
{qj}. For example, pairing between u and d quarks 
occurs for wd-Cooper pairs with any momentum in 
the set {q3}. To shorten notation, we will often refer 
to the A/ collectively as A and will henceforth write 
qf G {q/} simply as q". 

The free energy, f2, and the gap parameters. A, 
have been evaluated for many crystalline phases with 
condensates of the form for varied crystal struc- 
tures (i.e. varied {q/}) within the weak coupling and 
Ginzburg-Landau approximations A ^ (5/i <C /i [20| . 

In BCS-paired phases (like the CFL phase for ex- 
ample) quarks that pair have equal Fermi momenta 
even if in the absence of pairing their Fermi momenta 
would be split. This means that in such phases, the 
chemical potentials He, fJ-s and fis take on different 
values than in the absence of pairing, with the differ- 
ences being of order M^/fi. This rearrangement of 
Fermi surfaces exacts a free energy cost. The reason 
why crystalline phases can be more favorable is that 
they need pay no such cost: the quarks pair with- 
out rearranging Fermi surfaces; the gauge chemical 
potentials are as they would be in the absence of 
pairing, up to corrections that are proportional to 
A and thus negligible in the Ginzburg-Landau ap- 
proximation. We therefore set iie = M'^/{4:fi) and 
M3 = Ms = throughout, as in Ref. [20 ]. 

The free energy of a crystalline superconductor in 
the Ginzburg-Landau approximation relative to that 
for unpaired quark matter (i.e. the condensation en- 
ergy) can be written as 



2^2 

f^crystallino({ A/}) 



Y^Piai A}Aj 
. I 

\ / I>J ) 

- 5^7/(AJA,)3 + E 

77JJ AJA/A}AjA}Aj -|- 7123 A^Ai A2A2A3A3 



■0(A«) 
(13) 



r 



where Pi is the number of momentum vectors in been calculated for several crystal structures in an 
{q^}. The coefficients appearing in Eq. (I13|) have NJL model f29|, building upon previous analysis of 
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the two-flavor case [30|. The quadratic coefficient, turn vectors in the set {q^} and the Fermi surface 
a J depends only upon the magnitude of the momen- spUtting, S/j^j , and is given by 

I 



To find the most favorable magnitude of the vectors 
we need to minimize the free energy with respect 
to |qj| which at leading order is equivalent to min- 
imizing aj [30| . This fixes the ratio between the 
magnitude of the momentum vectors qf and to 
be the same for all the vectors, 

= |q?|=r/(5/i/ V/, a, (15) 

where 77 satisfies the relation, 

and is given numerically by 77 = 1.1997. 

The higher order coefficients in Eq. (fT5|) depend 
upon the relative directions {q/} of the vectors, 
and the task of finding favorable structures involves 
choosing polyhedra whose vertices are given by {q/}, 
calculating rJciystaiiine for each choice, and finding 
those which have the lowest ficrystaiuno- For a given 
choice of crystal structure, the quartic coefficients 
in ([in)) (i.e. the /3's) are dimensionless numbers 
times 5iJ? and the sextic coeflticients (i.e. the 7's) 
are dimensionless numbers times 5^i^. This makes it 
clear that the control parameter for the Ginzburg- 
Landau approximation is /^"^ /5^P = A^/(Mg /8/Lt)^, 
which is of order 1/4 for the most favorable crystal 
structures Because we shall only evaluate the 

phonon effective action and the shear modulus to or- 
der A^, we shall not need the values of the /3's and 
7's calculated in Ref. their role in the present 
calculation is only indirect, in the sense that they 
determine that the CubcX and 2Cubc45z structures 
arc most favorable. 

The free energy (|13p simplifies upon taking Ai = 
0. Ai can be neglected because it describes pairing 
between s and d quarks, whose Fermi momenta are 
separated by (5/ii which is twice 8^2 and (5 /is. 

The analysis of Ref. yielded qualitative ar- 
guments that two crystal structures called CubeX 



I 

and 2Cube45z should be most favorable, and ex- 
plicit calculation of their free energy to order A^ in 
the Ginzburg-Landau expansion showed that they 
have comparable free energy to each other and have 
a lower (i.e. more favorable) free energy than all 
the other crystal structures considered. Further- 
more, these two crystalline phases are favored with 
respect to unpaired quark matter and spatially uni- 
form paired phases including the CFL and gCFL 
phases over the wide range of parameters H]). For 
example, taking Aq — 25MeV and considering a 
fixed value of the strange quark mass Mg = 250McV, 
it turns out that one or other of the CubeX and 
2Cube45z phases is favored for ^ between 240 MeV 
and 847 MeV, which more than covers the entire 
range of densities relevant for neutron star interi- 
ors. We shall now describe the CubeX and 2Cube45z 
crystal structures. 



C. The CubeX and 2Cube45z structures 

The CubeX crystal structure is specified by two 
sets of unit vectors, {42} and {qg} depicted in the 
left panel of Fig. 1. Taken together, the two sets of 
vectors point towards the eight vertices of a cube. 
The four vectors in {q2} all lie in a plane and point 
towards the vertices of a diagonal rectangle of the 
cube, while the four vectors {qg} form the comple- 
mentary rectangle. We will use a coordinate system 
such that {q2} is given by {(1/V3)(±V2, 0, ±1)} 
(the four combinations of ± giving the four vectors 
in {q2}) and {qg} is given by {(1/V3)(0, ±^2, ±1)}. 

The 2Cube45z crystal structure is specified by 
two sets of unit vectors, {€[2} and {qg} depicted 
in the right panel of Fig. 1. The two sets {q2} 
and {qg} each contains eight vectors that point to- 
wards the vertices of a cube. The cubes speci- 
fied by {€[2} and {qg} are rotated relative to each 
other by an angle 45° about one of their C4 sym- 
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FIG. 1: (Color online) Left panel: The momentum vectors forming the CubeX crystal structure. This structure 
consists of eight vectors that belong to two sets {0(2} and {qs} which are shown as vectors which start from the 
origin. The four vectors in {q2} are given by {l/\/3){±\/2, 0, ±1) and point toward the vertices of the light shaded 
rectangle (pink online) that lies in the x — z plane. The four vectors in {qg} are given by (l/\/3)(0, ±\/2, ±1) and point 
toward the vertices of the dark shaded rectangle (purple online) that lies in the y — z plane. Taken together the two 
sets of vectors point towards the eight vertices of the light gray colored cube (only the edges are shown as light gray 
segments) . Right panel: The end points of the vectors forming the 2Cube45z crystal structure. This structure consists 
of sixteen vectors that belong to two sets {42} and {qa}. The eight elements of the set {q2} point towards the vertices 
of the black cube (only the edges are shown), and are given by (l/v'3)(±l, ±1, ±1). The eight elements of the set {qa} 
point towards the vertices of the light gray cube, and are given by {(l/\/3)(±\/2, 0, ±1)} U {(l/\/3)(0, ±\/2, ±1)}. 
The three dots denote the points where the axes meet the light gray cube, to clarify the orientation of the axes. 



mctry axis, passing through their common cen- 
ter. We will orient the coordinate axes such that 
{qj is given by {(1/V3)(±1, ±1, ±1)} and {q^} 
by {(1/V3)(±V2,0,±1)} U {(1/V3)(0,±V2,±1)}, 
which corresponds to a relative rotation by 45° 
about the z axis. 



The lattice spacing for 
crystal structure is [13, [lOl 



the face-centered cubic 



V3tt 4.536 36.28fi 



(17) 



For example, with M^/fi =100, 150, 200 MeV the 
lattice spacing is a =72, 48, 36 fm. The spacing 
between nodal planes is a/2. 



III. THE PHONON EFFECTIVE ACTION 



In this Section, we present our calculation of the 
effective action for the phonons in crystalline color 
superconducting phases of quark matter. In Subsec- 
tion A we describe the NJL model which we use. In 
Subsections B and C we introduce the phonon field 
and integrate the fermions out, yielding a formal ex- 
pression for the phonon effective action. In Subsec- 
tions D and E we introduce the Ginzburg-Landau 
approximation and evaluate the phonon effective ac- 
tion to order A^. 
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A. NJL model and mean field approximation 

As described in Ref. we work in an NJL 

model in which the quarks interact via a pointhke 
four-fermion interaction with the quantum numbers 
of single gluon exchange, in the mean field approx- 
imation. By this we mean that we add the local 
interaction term 



-'interaction 



■Gii^Tij)i^Tij) (18) 



to the Lagrangian density and treat it in the 
mean field approximation as we describe below. 
Here, we have suppressed the color, flavor and Dirac 
indices; the full expression for the vertex factor is 
T^^^^^ = Y{T^)ai35,j, where the are the color 
Gell-Mann matrices. To regulate the ultraviolet be- 
havior of the loop integrals, we introduce a cutoff 
A which restricts the momentum integrals to a shell 
around the Fermi sphere: /i^A< \p\ </i + A. 

The free energy Q, and the gap parameter A for 
a crystalline superconductor in the weak coupling 
and Ginzburg-Landau approximations depend upon 
the NJL coupling G and the cutoff A only via Aq, 
the gap parameter for the CFL phase with no Fermi 
surface splitting, given by 



Ao = 2 3 A exp 



(19) 



Model calculations suggest that Aq should lie be- 



tween lOMeV and 100 McV In the weak- 

coupling approximation, Aq ^ ^ and it is possible 
to choose A such that A ^ Aq while at the same 
time A <C /i. 

We will see that the explicit dependence on G dis- 
appears completely in our calculations and our re- 
sults depend on the overall strength of the interac- 
tions only via the values of the gap parameters A/, 
which in turn depend on G and A only via the combi- 
nation Aq [iOl . It will nevertheless be useful to keep 
the cutoff A in our loop integrals to sec explicitly 
that all the pieces involved are finite. The final ex- 
pressions for the coefficients of the terms that appear 
in the phonon effective action are ultraviolet safe, 
meaning, they are finite and A-independent as one 
takes A ^ oo. The weak-coupling approximation 
Aq ^ /i is crucial to this result, as this allows one 
to take A Aq while keeping A <C /it. Physically, 
the weak-coupling assumption means that the pair- 
ing is dominated by modes near the Fermi surfaces 
which means that the gaps, condensation energies 
and phonon effective action for the crystalline phases 
are independent of the cutoff in the NJL model as 
long as they are expressed in terms of the CFL gap 
Aq: if the cutoff is changed with the NJL coupling 
adjusted so that Aq stays fixed, the properties of the 
crystalline phases also stay fixed. 

In the mean field approximation, the interaction 
Lagrangian takes on the form 



J 



C 



interaction 



-3HvW + oV'^A(r)V' - -GiY{T^{V4')T){i'V), 



(20) 



where, tr represents the trace over color, flavor and 
Dirac indices, and where A(a;) is related to the di- 
quark condensate by the relations. 



A(x) 
A(x) 



7''At(x)7° . 



(21) 



B. Introduction of the phonon field 



We now consider the space- and time-dependent 
vibrations of the condensate, which will lead us to 
the effective Lagrangian for the phonon fields in the 
presence of a background condensate of the form 
(fT2|) . More precisely, we consider the condensate 



A(r) = AcF(r)®C7^ 



(22) 
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with 



in terms of which 



(r) 



7=1 



£lat3£lij^I 



q? 



Sometimes we will write 



with 



A/(r) = A 



7=1 



q? 



q/ 



(23) 



(24) 



(25) 



Note that / = 1,2,3 correspond to {ds), (us) and 
{ud) condensates, respectively. When we evaluate 
the phonon effective Lagrangian for the CubeX and 
2Cube45z crystals explicitly in Section IV, we will 
set Ai = and / will then run over 2 and 3 only. The 
condensate (j23p breaks translation invariance spon- 
taneously and we therefore expect Goldstone bosons 
corresponding to the broken symmetries, namely 
phonons. Phonons are small position and time de- 
pendent displacements of the condensate and, since 
the three condensates in can oscillate indepen- 
dently, we expect there to be three sets of displace- 
ment fields u/(x'). In the presence of the phonons, 
then. 



A,(r)->A^(x)=Az(r-u,(x)) 



(26) 



and we will denote the corresponding quantities ap- 
pearing on the left-hand sides of and as 
A(jp{x) and A"(x) respectively, i.e. 



Y.^io.pen.A'iix) , (27) 



7=1 



and 



(28) 



(We apologize that we have denoted the displace- 
ment fields, and hence quantities like A", by the let- 
ter u which in other contexts, but not here, denotes 
up quarks.) In the mean field approximation, the 
full Lagrangian, given by the sum of (|4]) and (|20|) . 
is quadratic in the fermion fields and can be written 
very simply upon introducing the two component 
Nambu Gorkov spinor 



X 



and hence x = ( "0 "0"^ 



(29) 



L = 



A-{x) (z^-^)^ j^ + 16G*'^^^ > ^ > 

(30) 

Here, = /i7o and /i is the matrix we have defined 
in Eq. pT|) . In the next Sections, we shall also often 
use the notation ffn = /.Ji7o, with i ~ 1,2,3 cor- 
responding to u, d, s respectively. The last term in 
Eq. (PD)) comes from the last term in Eq. which 
simplifies to (l/(16G'))tr((A")-'"A") for condensates 
given by Eqs. (P7)l and 



C. Integration over the x fields 

The spacing between vortices in the vortex array 
in a rotating superfluid neutron star is many mi- 
crons, and we will be interested in shear stresses ex- 
erted over lengths of order or longer than this length 
scale. This means that we need the effective action 
for phonon excitations with wavelengths of this order 
or longer. This length scale is many many orders of 
magnitude longer than the microscopic length scales 
that characterize the crystalline phase. For example, 
the lattice spacing is many tens of fm. This means 
that we need the effective action for phonons whose 
wavelength can be treated as infinite and whose en- 
ergy can be treated as zero, certainly many many 
orders of magnitude smaller than A. 

The low energy quasiparticles in a crystalline 
color superconductor include the displacement fields 
uj (x) , which are massless because they are the Gold- 
stone bosons of the broken translational symmetry. 
In addition, crystalline superconductors feature gap- 
less fcrmionic modes, as we now explain. In the ab- 
sence of pairing, quarks living at the Fermi surfaces 
can be excited without any cost in free energy; pair- 
ing in the crystalline phases yields gaps for quarks 
living in various ring-shaped bands around the Fermi 
surfaces, but leaves gapless fermionic modes at the 
boundaries of the pairing regions (loosely speak- 
ing, the remainder of the original Fermi surfaces 
other than the ring-shaped bands) [l^, [2^, [13, [3^ - 
The low energy effective theory includes fermions 
in the vicinity of the surfaces in momentum space 
that bound the pairing regions, in addition to the 
phonons that are our primary interest in this paper. 

To find the low energy effective action which de- 
scribes the phonons and the gapless fermionic exci- 
tations we need to integrate out those fermion fields 
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which have an energy larger than some infrared cut- 
off Air. For the application of interest to us, Air 
should be the energy of phonons with micron wave- 
lengths. If we were interested in thermal properties, 
Air would be of order the temperature T. (Either 
of these energy scales is ^ A, and by the end of this 
Subsection we will see that it is safe to set Air = 0.) 
In order to formally implement this procedure, we 
define 

ip = ip'j + tjji^ and hence, = ijj-^ + ^ , (31) 

where '(/'( and -0^ contain modes with energy in 
[0, Air] and ip-^ and -0) those with energy in [Air, oo]. 
Note that the boundary in momentum space be- 
tween the ij}) and modes will be nontrivial sur- 
faces that follow the boundaries of the pairing re- 
gions, where the gapless fermions are found. The 



corresponding decomposition for the Nambu Gorkov 
fields is, 

X = X)+ X{ and hence X = X) + X{ , (32) 

where X), X{i X) and X( are defined analogously to 
Eq. (pg)) . Carrying out the functional integral over 
the X) and X) fields will leave us with a low-energy 
effective action in terms of the u/, X{ and fields. 

We begin with the path integral expression for the 
partition function, 



(33) 



Z[n,Xi,X(]= J V[x)]V[x)Vf'''^^ , 



where the action of the Lagrangian in Eq. ((30)l can 
be written in terms of the decomposed fields (j32p . 
as follows, 



J 



16G 



tr((A")^A^' 



X) 



X) +X{S ^X{ 



(34) 



where the cross terms X)S and X{S "X) 
appear because the integration over imposes en- 
ergy conservation, and X) and X{ have support over 
different ranges of energy. The full inverse propaga- 
tor, S~^, in Eq. ((M]) is given by 



Since the Lagrangian is quadratic in the x and x 
fields, the standard result for fermionic functional 
integrals gives 



_ / i^ + lji A"(a;) 
" ' A"(a;) (i^-i^f 



(35) 



iS[u,X{,X{] = log(^[u,X(,X(]) 



XiS-'xi + ^tr((A^)^A«) 



-Tr„glog(5-i) 



(36) 



r 



where 5[u,X(,X(] is the low energy effective ac- 
tion that we are after, at present still given at 
a rather formal level. Here, Trng symbolizes the 
trace over the Nambu-Gorkov index, the trace over 
color, flavor, Dirac indices and the trace over a set 
of functions on space-time, with energies lying in 
[-oo,-Air]U [Air,oo' 



The factor i appears before 



Trjig because all the components of x and x are not 



independent. The actual independent fields are 
and ip. As promised, the effective action is a func- 
tion of the low energy quark fields, which appear in 
X(S'~^X(i and the phonon fields, which appear im- 
plicitly via the dependence of and A" on u/. 

We now concentrate on small displacements and 
hence drop all terms in the effective action of order 
(u/)'^ or higher. This is most simply done by looking 
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at 



1=1 

3 



q" 



(37) 



r 



where we have defined 



2q?.U7(x) =0?(x) 



(38) 



We will refer to both the uj fields and 0° fields as 
phonons, as we can write one in terms of the other. 

We now argue that as far as the calculation of 
the shear modulus is concerned, we can look only 
at the part of the effective action that describes the 
phonons. The remainder of the effective Lagrangian, 
where the low energy quark fields appear, can be 
written as 



with 



and 



X{S = ^f+^f<p 
i^ + l/i A(r) 



(39) 



^/=X((^l(rf [^f-'^^)T)xi (40) 



r - ( A"(x)-A(r) , 

(41) 

We shall sec in Section HV Al that the shear modulus 
is related to the coefficient of 



(42) 



dx"" 



in the Lagrangian, which makes it obvious that £/ 
does not contribute. The coefficient of (|^^ at the 
scale Air receives contributions from the X) ^-nd 
X) fermions which have been integrated out. The 



phonons and low energy quarks (at an energy scale 
lower than Air) interact via the Cf^ term in the 
Lagrangian. Formally, then, one has to solve con- 
sistently for the phonon propagator and the quark 
propagator, which are coupled. However, the ef- 
fect of the phonon-fermion interactions on the self 
consistent calculation of the gap parameter using 
£/ is small, because the quark loops come with 
an additional factor of /i^, which is is much larger 
than (Air)^, and hence the quark propagator can 
be considered to be unaffected by the phonons. 
The phonon propagator, and hence the shear mod- 
ulus, will depend on the phonon-fermion interac- 
tions, meaning that the phonon propagator and con- 
sequently the coefficient of (|42|) will run as Air is re- 
duced. However, as long as Air is much smaller than 
A, |q| and the change in the value of the shear 
modulus from integrating out more fermions below 
the scale Air will be negligible compared to its value 
at Air . This means that we can take Air = 0, inte- 
grating all of the fermions out from the system and 
obtaining an effective action for the phonons alone. 
This procedure is correct for the calculation of the 
shear modulus but would not be correct for, say, 
calculating the specific heat of the system, which 
is dominated by the gapless fermions not by the 
phonons. We have checked numerically that the dif- 
ference between the shear modulus calculated with 
Air = and that with a small but nonzero Air is 
negligible. 

Finally, therefore, the effective action we are in- 
terested in depends only on the phonon fields, and 
is given by 



i5[u] = log(Z[u]) 



16G 



tr((A")'^A") 



-Trnglog (5- 



(43) 
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where now the Tr,ig includes a trace over functions 
in space-time containing all energy modes. 

For the single plane wave "crystal" structure in 
which only one of the A/ is nonzero and {q/} con- 
tains only a single wave vector [2^, [2^ , we can in- 
vert the Nambu-Gorkov inverse propagator in the 
absence of phonons without expanding in A, and can 
therefore obtain the effective action for the phonons 
up to second order in (j), to all orders in A. We shall 
do this exercise in Appendix |B] For the realistic 
crystal structures, CubeX and 2Cube45z, however, 
we cannot invert the full inverse propagator and we 
therefore proceed by making a Ginzburg-Landau ex- 
pansion in A. 

D. Ginzburg-Landau expansion 

In order to obtain the effective action for the 
phonon field we first separate the full inverse prop- 

I 



where, the Tr on the right hand side is over Dirac, 
color, flavor and space-time, and we have used the 
cyclic property of the trace to equate the two con- 
tributions obtained from the trace over the Nambu 
Gorkov index. Finally, substituting (1^5]) back in P5)) 



agator, S~^, defined in Eq. ((35l) . into the free part, 
5*0^^ and a part containing the condensate, S, as fol- 
lows: S^^ = Sq^ + S, where 



and 



Then, we can expand the term log(S' ^) that ap- 
pears on the right-hand side of Eq. (|43p as 



(46) 



(47) 



(48) 



(49) 

I 

and simplifying the Dirac structure of the operators 
using the Cj^ Dirac structure of the condensate and 
the properties of the charge conjugation matrix C, 
namely C(7^)-^C~-^ = —7'' and = —1, we obtain 



Tr„g(log(5o-i + S)) = TvneilogSo') + Tr„g(S'oS) - iTr„g(S'oS)' + ... 
where Trng(logS'(^^) is related to the free energy of unpaired (Normal) quark matter fl^ hy 

Trng(log%i) = -i2nN J d'^x = -i2VTnN , 
with VT the space-time volume. Since 

(^oS) = 

only even powers of contribute to the trace over Nambu Gorkov indices and we can write, 

Tr„g(log(5-i)) = -i20w(FT) - -Tr((*^ + i^Y^ /^^{x){{i^ - i^f^A^i^ 



_ (i^ + /i)-iA"(x) 
((i^-/i)^)-iA"(x) 



n 

n=l 
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S[u] 



AG 



11 

n^iVT) -Tr((*^ + i^)-'A"cFix)iii? - 



2i ' — ' n 

n=l 



where the trace tr cf is over color and flavor indices 
and where A'^p{x) depends on u/(a;) via Eqs. ([57]) 
and ([55)1 . Eq. ([5D|) is the low energy eff'ective ac- 
tion for the phonons. written as a Ginzburg-Landau 
expansion in A. We will calculate the leading con- 
tribution to S[u], namely that proportional to A^. 

The first term on the right hand side of Eq. ([5(1)) 
does not have any derivatives acting on u/ and hence 
can only contribute to the mass of the phonon. which 
we know must be zero by Goldstone's theorem. In 
Appendix |^ we show explicitly that the uj depen- 
dence in / d'^xtTcF{{AcFV Aqp) cancels out, and 
its value is given simply by 

^ / d'xtrcF{{A"cF)^AhF) 

1 (51) 
^{VT)-Y,{AjA})Pj , 

I 

where Pj is the number of plane waves in {q/}. 

We now proceed to evaluate the third term on the 
right hand side of Eq. ([5(1)1 diagrammatically. We 
will expand the action S[u] in Eq. ()50p in powers of 
(j) (or equivalcntly U/) up to second order in (f> by 
using the Feynman rules described in Fig. |5| 

The lowest order term is independent of (j). The 
sum of this term and of the term given in Eq. ([5ip 
has a simple interpretation. In the absence of 
phonons, and considering that the fermionic fields 
have been integrated out, the action in Eq. ([50)1 



turns out to be proportional to the free energy of 
the system. More specifically, 

S[U = 0] = -(1/T)(r!erystallinc + ^n) , (52) 

where rJcrvstaiUno is given as a Ginzburg-Landau se- 
ries in A [201 • Since ([5^ is independent of (j) it does 
not affect the equations of motion of the phonons 
and we will simply drop it from our calculation. 

We now consider the term that is linear in (p. We 
will evaluate the leading term in the action propor- 
tional to A^, which we will call S'^^ and which is 
represented diagrammatically in Fig. |3) Both terms 
shown in Fig. O are proportional to the trace of 
£ia/3£ja/3£iij£jij, which is nonzero only if / = J 
and therefore only terms proportional to AjAj are 
present. We could have anticipated this result from 
the symmetries of the problem. The Lagrangian 
conserves particle number for every flavor of quarks, 
which corresponds to symmetry under independent 
global phase rotations of quark fields of the three 
flavors, meaning independent phase rotations of the 
three A/. The effective action should be invariant 
under these rotations and hence A/ can only occur in 
the combination AjAj. (Although the condensate 
spontaneously breaks them, the requirement is that 
the Lagrangian has these symmetries.) Then, the 
sum of the diagrams in Fig. [31 which corresponds to 
the contribution to the action linear in the phonon 
field and second order in A, is given by 



S^-' =Y^A}A,J: EJ^.J (g,(2.)V^n2q?-2qJ + .) 



tr 



0?(fc)-0^(fc) 



(53) 



I 



where k is the four momentum of the phonon field on the right-hand side ensures momentum conserva- 
and the trace is over Dirac indices. The Dirac delta 
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FIG. 2: Propagators and interaction vertices for the Lagrangian up to order </i^. The dashed lines represent propa- 
gating quarks, and the solid lines represent propagating quark holes. The dot-dash lines represent external phonons. 
The subscript on is the index of the flavor which is propagating, and determines the value of the chemical potential 
that appears in the propagator. The Aj vertex comes along with a momentum insertion 2qf and a vertex factor 
£iai3£iij- Similarly, Aj comes with a momentum insertion — 2qJ and the same vertex factor. 



tion, 



2qJ - 2q| + A: = 



(54) 



meaning that the net momentum added to the loop 
is zero. But since we are looking at the low energy 
effective theory, we can take k much smaller than 
the momentum vectors qj and therefore Eq. ([5^ 
can be satisfied only if fc = and qJ = qj, which 



means that a 
vanishes: 



h. Using this resuh, we find that 



S"^^' = 



(55) 



That is, the term linear in cj) is absent to order A^. 
One can similarly argue that it is absent to all orders 
in A. 
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i>H-m-i,)-ii^k 



FIG. 3: Diagrams that contribute to order (^A^. The dashed lines represent propagating quarks, and the sohd 
hnes represent propagating quark holes. The interaction vertices have been defined in Fig. [5] and the color-flavor 
structure is also indicated. Note that the trace over the color-flavor epsilon tensors, eiapejafsenjeiij forces I = J 
and momentum conservation implies qf = qj, as well as fc = 0. The two contributions are then equal in magnitude 
and opposite in sign, and hence cancel. 




FIG. 4: Diagrams that contribute to order (fi^A'^. In drawing the diagrams, we have used the fact that the trace of 
the color-flavor tensor forces I — J. We have also used the fact that momentum conservation requires that the net 
momenta added by the condensate and the phonons are separately zero. In the first diagram, momentum conservation 
at the A| vertex imposes ^2 — fci = k. 



Now we consider the terms of order 0^ , which give 
the first nontrivial contribution to the phonon effec- 
tive action. We will evaluate these terms to order 
and we will call the corresponding contribution to 



the action S'^ . The terms contributing to S'^ 
arise from the diagrams given in Fig. [4l and give 



EE 



(2^)4 



(56) 
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where fc = fc2 — fci is the four momentum of the phonon and 



tr 



(57) 



where the trace is over Dirac indices. In the next 
Subsection we evaluate Eq. ([57]) . The reader not 
interested in the details of our calculation will find in 
Eq. ([75]) the final expression for the effective action 

E. Evaluation of 5'*''^' 

We turn now to the evaluation of Vf{k) of Eq. ((57)) 
and hence, via (|56p . the loading nontrivial contri- 
bution to the phonon low energy effective action, 

To begin, if we set fci = fc2 = in ([57]) . implying 
that fc = 0, we see that 7^/(0) = 0. In this way, we 
see explicitly that the phonons are massless to order 
A^. As mentioned before, we are interested in the 
low energy, long wavelength phonons. We therefore 
expand Pj{k) in powers of k and drop terms of order 
and higher. 

We are working in the limit in which Sfi, q = |q| = 
T]6fi and A are all much smaller than (A ^ /i 
follows from the weak coupling approximation and 
(5/i ^ /.t follows from requiring The 
Ginzburg-Landau approximation, which is the fur- 
ther requirement that A^ <|; 6^^, is not required in 
the derivation of the simplifications of Eq. (|57p that 
follow.) This means that we can choose the ultravi- 
olet cutoff A such that 6^, q, A <^ A <^ fi. Since the 
largest contribution to the integrals comes from the 



region close to the Fermi surfaces and since A ^ /i, 
the integration measure in Eq. (|57p can be approxi- 
mated as follows: 

where v is the unit velocity vector in the direction of 
the spatial momentum vector, v = p/|p|, and J dv 
represents the integral over solid angle covering the 
Fermi surface. The residual momentum s is defined 
by the relation s = |p| — /I, where ft is an energy scale 
lying close to the quark Fermi surfaces. In evaluat- 
ing Vf{k), we will take p, to be the arithmetic mean 
of and i^ij, for convenience, but since the two 
integrals on the right-hand side of Eq. (|57p go as 
log(A) for large momenta, choosing any other value 
for fL close to the Fermi surfaces changes the value 
of V^k) only by 0{S^/A), which we will ignore. 
We introduce the two null vectors, = (liv) and 
= (1, — v), as is done in the High Density Effec- 
tive Theory (HDET) It is also useful to define 
four momenta = (p'',sv) = p^ — (0,/iv), which 
can be thought of as residual momenta as measured 
from the Fermi surface. It is easy to verify that 
V ■ I = p'^ — s and V ■ I — p^ + s. In the weak cou- 
pling limit, for a generic four vector p' that is small 
compared to A, the propagators in Eq. (|57p simplify 
as follows: 



19 



and, similarly. 



(/ + (p')"+A*,)7"-(P + P')-7 



i> + i>' + i/ij (pO + + H - IP + P'I)(P" + b')" + N + |P + P'l) 



M7 - P ■ 7 



7 — V • 7 



2 ^vpO + {p'f - s + - /i) - p' • V 
1 / T/-7 



2 V^-(/+p') + (Mj-A). 

(/ + (pO"-Mfe)7°-(P + PO-7 

i> + i>'- h iP° + iP')° -f^k-\p + p'Db" + (P')° - Mfc + |P + P'l 

7" + V • 7 



1 



_ 1 

^ 2 yp" + {p'f + s- [fik - A) + p' ■ V 

2 \v.{l+p')-{^^-fi)^ 
Upon using these simplifications in Eq. (|57p . we obtain, 

r d»" r dv 1 



_^ 2mJ_f, J ^t: {V ■l~ir-cfj+5^ii){V -l-w-cfj+Spi) 



dp;; 

J-K y 47r • (/ + fci) - V • qf + • (/ + fca) - v • qf + 5pi) 



(dm 



(59) 



(60) 



(61) 



By making the changes of variables p'^ — > — 
s — > — s, V —V and J/i/ — > —dfii it is easy to 
show that 7'^(-fc°,k) = rf{k^,k). In addition, it is 
clear from Eq. ((56)) that only the part oi'Pf{k) that 

,2 A 2 

is even under fc — > — fc contributes to 5''' . Hence, 
terms in the expansion of T'f (fe) proportional to odd 
powers of k do not contribute to the effective action. 
Furthermore, at second order in k this implies that 
there cannot be terms proportional to k^kcf)^/^'^ in 
the effective action. (Terms like fcgk^^^A^ can of 
course appear, but are higher order in k.) This is 
useful because we can handle the spatial and time 
components of k independently, thereby simplifying 



the calculation of Vfik). 

In order to simplify the calculation we rewrite 
'Pj{k) a little differently. Multiplying and divid- 
ing the integrand appearing in the second term in 
Eq. (|6T|) (the integrand depending upon ki and /C2) 

by 

{v■{l + k,)-w■ci^, + 5^lI) 

X (F • (/ + fc2) - V • q? + 5^iI) ^ ' 

and collecting the term with numerator 

{V-k){V-k) = (fc2-fci))(l7. (fc2-fci)), 

after some algebra we obtain 



2 {V-k){V-k) 1 

{V ■{l + ki)-w-cfj+ 5ni){V • (/ + /C2) - V • qj + 5ixi) D{1 + ki)D{l + k2) D{1 + ki) Z)(/ + fc2) 



k2) 
(63) 
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where 



D{1) = {V-l + S^ii - V • q?)(F • / + S^iI 
We can then write "P/ (fc) as Vfik) ^ la + h with 



^ 2 7r2 / 27ri 



2 7r2 



dv {V ■k){V- k) 
_A J D{l + ki)D{l + k2) 
rfv 1 1 /i^ /■ (ip° 

2TTi 



27ri 



2 7r2 
1 



{5^.11 



ds 



dv 1 



4^ D(Z + fca) 
I 



_A J AT:D{l + ki) 



(64) 



(65) 



(66) 



The reason to separate Vf^k) as a sum of /q and /i 
wiU become clear in a moment, when we argue that 
/o = 0. 

We now proceed to evaluate /q and Ii separately. 
As we discussed, we can consider the spatial and the 



temporal parts of k independently, and we begin by 
taking k = (0,k) with k = k2 — ki. In this case Iq 
can be expressed as a sum of three terms, each pro- 
portional (up to a prefactor /tt^ or —ii^/{2Tr^)) 
to an integral of the form 



n(q?,5A^,,p')= / ^ 
27ri 



ds 



-A 
A 



dv 



1 



ds 



-A 



47r D{l+p') 
dv 



(67) 



47r (pO 



V • p' - V • + 5^i){p^ + ,s + V • p' - V • qj + (5/^/) 



where p' is in this case a purely spatial vector, p' = (0,p'). This integral can be evaluated by following the 
steps outlined in Ref. [s^l- We first perform a Wick rotation, p° ip'^, and then do the p^ integration by 
the method of residues, followed by the ds and dv integrals. For p' = 0, the integral is calculated in Ref. |30l | 
and is given by 



n(q?,<5M/,0) 



1 



log 



A2 



(q«)2-<5/i2y 2|q- 



— r log 



\<fI\-s^lI 



(68) 



By making the change of variables s — > s — v • p' in Eq. (j67p . we see that the integrand appearing in the 
definition of n(qj, (5/x/,p') in Eq. ([67)) . can be written as the integral appearing in n(q°, (5^/, 0), but with 
the limits of s integration changed to [—A — v • p', A — v • p']. Since, 11 goes as log(A), however, this change 
in limits changes the value of the integrand only by a quantity of order k/A, which we ignore. Thus, 



T — ^ 



n(q?, (5^/, 0) - ^n(q?, Sfij, 0) - in(q?, (5^,, 0) 



(69) 



The integral /i can be evaluated using similar steps, namely Wick rotate p^ — *■ ip"^, perform the p^ integration 
by residues and then do the ds and dv integrals. The final result is 



n(q?, ,5^^/, 0) - in(q? - ik, s^lI, o) - in(q? + ik, s^u, o) 



(70) 
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We note that the final result depends only upon k = k2 — ki . Expanding Ii in k, we find 



1 



4(qf 
1 



4(qJ)2 



1 



6fij) 2(q?)2 V 2|q? 
SfJ-i f\ q,j\ - 

<5^/ 



— r loe 

a o 



r 



\<fi\~Sfii 
|q?l + '^A^/ 



(71) 



where k||^ is the component of k which is parallel 
to qf , k||a = qj(k • qf ), and k±a is the component 
perpendicular to q", kj^^ = k — kn^. In deriving 
Eq. (j7ip we did not assume any particular relations 
between q" and Sfii, but now we choose the value 
of given by Eq. that minimizes the free 
energy. Substituting Eqs. US]) and ^ into Eq. ^ 
simplifies Vf{k) considerably, yielding 



Vfik) 



2 



1 



iSfijir]^ - 1) 



for k = (0, k) , 
(72) 

where we have dropped the terms of order fc^. 

The final expression for Vj{k) has the following 
features. First, the k_Lo term has dropped out. 
This means that for a single plane wave condensate, 
phonons that propagate in the direction orthogonal 
to the plane wave that forms the condensate cost no 



energy up to order (kj^a)^A-^0^. Second, we note 
that the coefficient in front of (k||a)^ is negative. 
This means that the crystal structure is stable with 
respect to small fluctuations in the direction of q°. 
(Recall that action goes like kinetic energy minus 
potential energy; since here k is spatial, we have 
potential energy only meaning that decreasing the 
action corresponds to increasing the energy, hence 
stability.) This result is a direct consequence of the 
fact that we chose jq^j to minimize a/, meaning that 
any deviation from the most favorable modulation of 
the condensate in the direction of qj increases the 
free energy of the system by an amount of order . 



We now evaluate /q and Ii in the case where k 
is purely temporal, namely k = (fc*',0) with fc" = 
fcj — fc]'. With these fci and ^2, the value of Ii turns 
out to be 



T ^ 



1 / (JM/-fcg + (fc72))^-(q?)^ A ^ . . _ ,0 ^ .,0/2^^ w / ^A^. - fe? + (^72) + |q°| 
2^°H (5/.. - fc?)2 - (q?)^ ^i + ^^/')^^°H^A^/-fc? + (fc°/2)-|q?i 



^ (6ur fc") lo. f^J^L^RlM] + 1 W A'?M/-fc2"-(fe72))^-(qzf \ 



,5Ai/-fcO + (fcO/2)-|q?|. 



;.0 / 7 n 



(fc"/2))log 



6ni - k^- {k"/2) + \q]\ ^ _ / Ofii - «.-2 -r iqj 



(5^lI - kl + \cfj\ 



(73) 



Although, it appears 
we find. 



that Eq. ([73]) does not depend solely on fcj ~ ^ii upon expanding in small fej and ^2, 



h = iky ^ 



^2 4((q?)2 - ^^2) 
Turning now to /q, this quantity is given by 



+ 0(fc4)« (fcO) 



,0\2 



7r2 4,5Ai?('f - 1) 



in — ^ 



n(q?, <5Ai7, 0) - ^n{q'}, s^iJ, k,) - in(q?, (5/.,, fc2) 



(74) 



(75) 
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where ki = {k^, 0). When its third argument is a purely temporal four- vector, 11 is given by 

1 



27ri 



-A 



47r + p'O - s - V • qf + 5^ii){p° + + s - w ■ cfj + Sfij) 
I 



(76) 



W) that 
/o 



where p' ~ {p"^ , 0). It is apparent from Eq. 
by making the change of variables p^ —> p'^ + p'" we 
obtain Il{c[j,d^i,p'^) = n(qf , 0), leading us to 
conclude from Eq. ([75]) that /q = 0. We advise the 
reader that in order to obtain the result /q = by 
the approach that we have employed in a straight- 
forward manner, it is important to shift po before 
Wick rotating. (Note that if we calculate Iq using 
dimensional regularization we find /q = in agree- 
ment with what we obtained by change of variables.) 



Finally, therefore, with /q = we obtain 



vnk)^ikr-. 



7r2 AS^ijiif - 1) 



for k = (A:",0) 



(77) 

We see that this comes with a positive sign, as is 
appropriate for a kinetic energy term. 

Substituting the expressions given in Eqs. (|72p 
and ([77)1 back into the action ([55)) . we obtain 



EE 



d'^k 



^^^y^kWAk) [kl - (k||,)2] ^ .^^^^^^^ _ 



(78) 



where k||a = q/(k • qj) and (/^/(fc) = u/ • (2qj). Inverse Fourier transforming back to position space, and 
taking out a factor of half for future convenience, we obtain the effective action for the displacement fields: 



5[u] 



/ ' q? 



u,)ao(q? • u,) - (q? • a)(q? • u,)(q? • 9)(q? • u,) 



(79) 



This is the low energy effective action for phonons 
in any crystalline color superconducting phase, valid 
to second order in derivatives, to second order in 
the gap parameters A/'s and to second order in the 
phonon fields ui . This is the central technical result 
of our paper. 

Because we are interested in long wavelength, 
small amplitude, phonon excitations, expanding to 
second order in derivatives and to second order in 
the ui is satisfactory in every respect. Not so for 
the expansion to second order in the A/. As we 
have discussed previously, the Ginzburg-Landau ap- 
proximation is at the point of breaking down in the 
most favorable CubeX and 2Cube45z crystal struc- 
tures. Before proceeding, we therefore ask what kind 
of corrections to ([75)1 will arise at higher order in A. 
The first thing to note is that in the weak coupling 
limit /i appears only as an overall factor of /i^ in 
front of the fermion loop integrals. After simplify- 



ing the fermionic propagators as in ([59|) and ([60]) 
and taking A to oo, the only two independent di- 
mensionful quantities that remain in the integrals 
are k and (5/i/. (Recall that |qj| is given by jyiJ/i/ 
and so is not independent.) Since we found the ac- 
tion only up to terms which are second order in the 
derivatives and second order in A, to ensure the La- 
grangian density has dimension four, only a dimen- 
sionlcss factor can multiply /x^jA/pS^Uj, as we can 
sec is true in Eq. (|79p . Higher powers of A^ will 
appear in Eq. (|78p in combination with higher com- 
pensating powers of 5pL~^ . Consequently, there will 
be corrections to the coefficients of fcg and {kwaY 
in ([75)1 suppressed by factors of (A^/J/i^) relative 
to the leading order result that we have obtained. 
In addition, there will be new terms. There is no 
reason to expect that the coefficient of (k^a)^ will 
remain zero at 0{iJ?\/!i^i\'^{du.iY / 5^j?). Finally, we 
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see that there are no terms in (|78p that "mix" the 
different u/(fc). This follows from the color-flavor 
structure of the condensate as discussed above. At 
higher order, there will be terms proportional to 
fi^\AiAj\^duiduj/Sn'^, which do "mix" the differ- 
ent u/'s. 

With the phonon effective action now in hand, 
in Section IV we shall relate the coefficients of the 
terms in S(u) involving spatial derivatives acting 
on the displacement fields to the shear modulus of 
crystalline color superconducting quark matter with 
specified crystal structures. 



IV. EXTRACTING THE SHEAR MODULUS 



We see from Eq. ([79| that the action of the phonon 
fields, 5(u), is a sum of two terms: the kinetic en- 
ergy, which has time derivatives acting on the fields 
u, and the potential energy, which has spatial deriva- 
tives acting on u. From the basic theory of elastic 
media [6l| , the potential energy is related to the elas- 
tic moduli that describe the energy cost of small de- 
formations of the crystal. In this Section, we present 
this relation explicitly and calculate the shear mod- 
ulus for the CubcX and 2Cube45z crystal structures. 



A. Generalities 

Let there be a set of displacement fields u/ prop- 
agating in a crystalline color superconducting ma- 
terial. (We will set the problem up in the general 
case where all the A/ are nonzero, meaning that / 
runs from 1 to 3.) The kinetic energy density for the 
displacement fields takes the form 



K. = 



(80) 



IJ run 



where u^* , u j are the space components of the vec- 
tors u/ and uj respectively. (We will use the in- 
dices m, n, u and v to represent spatial indices in 
the following). As we are working only to order A'^, 
the only nonzero components of are those with 
I = J. We will choose the direction of the axes, x, y 
and z such that for every / and J, p™/" is diagonal in 
the m and n indices and we will denote the diagonal 
components of p™" by p™. We can then rewrite the 



kinetic energy density as 



At higher orders in A^, we could need to choose 
a new linear combination of fields = AJ^J^u" to 
render the kinetic energy diagonal in the IJ and mn 
indices. 

The potential energy density to quadratic order in 
the displacement fields can be written as 



VflUW 

IJ 



IJ 



^ duf dvJ} 

dx^ dx"" 



(82) 



where A"™"" is the elastic modulus tensor 



The com- 
ponents of the tensor du™/dx^^ that are antisym- 
metric in the mn space indices are related to rigid 
rotations. The symmetric components of the tensor, 
namely the "strain tensor" 



duf duf 
, — — H 

2 V dx"" 9x™ 



(83) 



tell us about deformations of the medium. In the 
previous Section we have shown that, to order A^, 
there is no interaction between the displacement 
fields u/ and uj with / different from J. There- 
fore A^j""" is diagonal in the / and J indices and, 
denoting the diagonal entries by A™""", we find 



jEE^r 



(84) 



Next, we define the stress tensor acting on the 
crystal / as 61| 



dU 

dsf^ 



(85) 



which is symmetric in its spatial indices. For a po- 
tential U that is quadratic in the displacement fields 
and is given by the stress tensor is 

erf" = Xf^^^'sY . (86) 

The diagonal components of a are proportional to 
the compression exerted on the system and are 
therefore related to the bulk modulus of the crys- 
talline color superconducting quark matter. Since 
unpaired quark matter has a pressure ~ /i^, it gives 
a contribution to the bulk modulus that completely 
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overwhelms the contribution from the condensation 
into a crystaUine phase, which is of order /i^A^. We 
shall therefore not calculate the bulk modulus. On 
the other hand, the response to shear stress arises 
only because of the presence of the crystalline con- 
densate. 

The shear modulus can be defined as follows. 
Imagine exerting a static external stress ct/ having 
only an off-diagonal component, meaning cr™" 7^ 
for a pair of space directions m ^ u, and all the other 
components of a are zero. The system will respond 
with a strain s™ satisfying ([55]) . The shear modulus 
in the mu plane, , is defined as half the ratio of 
the stress to the strain: 



a quadratic potential, with cr™ 
shear modulus is 



given by ([86)) . the 



\raunv „nv 



(88) 



where n and v are summed but m and u are not. 
For all the crystal structures that we shall consider 
below, the only nonzero entries in A™""^ with m ^ u 
are the A™"™" entries, meaning that 186]) simplifies 
even further to 



1a 



mumu 



(89) 



257^" 



(87) 



where the indices m and u are not summed. For 



again with m and u not summed. 

Putting Eq. ([H]) and Eq. ([84|) together, the action 
for the displacement fields can be written as 



(90) 



The equations of motion obtained by extremizing 
the action S[u] with respect to the displacement 
fields u are 



Pi 



(91) 



where I and m arc not summed. The dispersion 
relations are found by solving 

Det [pTklSmn - XT'^'kuK] = (92) 
for all /, where m is again not summed. 



B. Elastic moduli of crystalline phases 



In order to set up the extraction of the elastic 
moduli of crystalline phases, we need to rewrite the 
action ([79]) for a generic crystalline phase in a form 
which makes comparison to (j90p straightforward. 
Writing the spatial indices in Eq. ([75)1 explicitly, we 
obtain the low energy phonon effective action in the 
form 



J 



5[u] = i J d^xJ2^j 



(gj)'"(g?)" I (aour)(9ou7) - I 1 {duuTm.n^) 



r 



(93) 



where we have defined 

'^^ ~ 7r2(?72 _ 1) 



For a given crystal structure, upon evaluating the 
sums in and then using the definition ([55)) to 



(94) 
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compare (|93|) to ([90]) , we can extract expressions for 
the A tensor and thence for the shear nioduh. The 
quantity k is related to the elastic modulus for a 
condensate whose "crystal" structure is just a single 
plane wave, as discussed in Appendix IB] 

In the next two Subsections, we will calculate the 
shear modulus for the CubeX and 2Cube45z crys- 
tals. We will not discuss the expression for the ki- 
netic energy density, /C, but it is easy (and necessary) 
to check that in each case below we have chosen our 
axes such that /C only contains terms that are diag- 
onal in the spatial indices m and n. 

Note that henceforth we set Ai = and A2 = 
A3 = A, meaning that 



^2 



K-3 



2^^|Apr/2 
7r2(?/2 — 1) 



0.664^^1^2 



(95) 



The shear moduli that we evaluate each take the 
form of a dimensionless constant times k. 



C. Shear modulus for the CubeX crystal 



Orienting the axes as shown in the left panel of 
Fig. 1, we have {qa} = {(1/V3)(±a/2, 0, ±1)} and 
{qg} = {(l/y3)(0,±V2,±l)}. Calculating the rele- 
vant sums and substituting in (|93p . we find that the 
potential energy is given by 




Recall that the only components of the stress ten- 
sor that arc relevant to the calculation of the shear 
modulus are given by dU/ds™^ for m 7^ u. These 
are 
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'2 ^ "2 — "^'^ ^2 



and 
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'3 — ^3 



from which we obtain 



'3 ' 



and 



2s^^ 9 



yz zy ^3° 16 

3 3 2s.y 9 



(97) 



(98) 



(99) 



(100) 



We can display the result succinctly by writing 
two shear matrices V2 and 1/3, which have only off- 
diagonal entries and are symmetric in the spatial 



indices: 



1^2 




V3 




(101) 

The zeroes in these matrices are easily understood. 
The A2 crystal is translation invariant in the y- 
dircction, because all the wave vectors in the set 
{q2} lie in the xz-plane. This means that the xy- 
and yz-components of V2 are zero. The only nonzero 
shear modulus is that for shear in the xz-plane. Note 
also that the A2 crystal has nonzero A|^^^ and A|^^^ , 
meaning that it has a nonzero Young's modulus for 
compression or stretching in the x- and z-directions 
confirming that, as the shear modulus indicates, it 
is rigid against deformations in the xz-plane. Sim- 
ilarly, the A3 crystal is translation invariant in the 
x-direction, meaning that the only nonzero compo- 
nent of the shear modulus is that for shear in the 
yz-plane. 

The vortices in rotating crystalline color super- 
conducting quark matter have currents of m, c? and 
s quark-number flowing around them, meaning that 
the phase of both the A2 and A3 condensates winds 
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once by 2tt around a rotational vortex, and meaning 
that both A2 and A3 vanish at the core of the vor- 
tex. This in turn means that it will be free energeti- 
cally favorable for the vortices to be pinned at places 
where the A2 and A3 condensates already vanish in 
the absence of a vortex. The A2 crystal has two 
families of nodal planes where A2(r) vanishes. One 
class of nodal planes are parallel to the ly-plane and 
are located at z = [(2n+l)'K^/?>) / {Aq), where n is an 
integer. The others are parallel to the yz-plane and 
are located a.ix = ((2n+l)7rV6)/(4q). Similarly, the 
A3 crystal has nodes along z = ((2n -|- l)7rV3)/(4(j') 
and y = {{2n + l)7rV6)/(4g). So, we expect that 
the most favorable location of the vortices will be 
within the common nodal planes of the A2 and A3 
condensates, namely, z = ((2n -I- l)7rV3)/(4g). If 
these vortices are oriented in the x-direction, they 
will preferentially be located (i.e. will be pinned at) 
at x = ((2n + l)7rV6)/(4g). And, if the vortices in 
an array of vortices oriented in the x-direction try to 
move apart (i.e. move in the yz-plane) as the rota- 
tion slows, in order to move they will have to shear 
the A2 crystal which has a nonzero I'l^. Similarly, 
if the vortices are oriented in the y-direction, they 
will preferentially be located (i.e. will be pinned at) 
at y = ((2n -I- l)Tr^/6)/ {4q). And, if the vortices in 
an array of vortices oriented in the y-direction try 
to move apart (i.e. move in the xz-planc). they will 
have to shear the A3 crystal which has a nonzero 
j/f^. Thus, the nonzero shear moduli that wc have 
found in (|10ip arc sufficient to ensure that vortices 
pinned within the CubeX phase are pinned to a rigid 
structure, with the relevant shear modulus having a 
magnitude 16k/9. 

We note as an aside that further evidence for the 
rigidity of the CubeX crystal can be found by eval- 
uating the phonon velocities and showing that at 
long wavelengths the velocity of transverse phonons 



(which are found in a rigid solid but not in a fluid) 
is comparable to that of the longitudinal phonons 
which are found in both fluids and solids. We will 
evaluate the velocities of the longitudinal phonons 
upon ignoring the existence of longitudinal oscilla- 
tions in the gapless fermions, which have velocity 
1/a/3 in the limit of weak coupling. For this reason, 
the longitudinal phonon velocity that we calculate 
should be seen only as a benchmark against which 
to compare the transverse phonon velocity. The true 
sound modes would be linear combinations of the 
longitudinal phonons and the fermionic sound waves, 
which must in reality be coupled. This complica- 
tion does not arise for transverse phonons: the fluid 
of gapless fermions has no transverse sound waves; 
they can only arise as excitations of a rigid structure, 
like the crystalline condensate wc analyze. Consider 
as an example the phonons of the A2-crystal. From 
the dispersion relations ([5^ it is easy to show using 
pf = 8k/3, P2 = and P2 = 4k/3 that longitu- 
dinal phonons propagating in the x-direction have 
V ~ ■\/2/3 while transverse phonons propagating in 
this direction have the same v = y/2/3. For prop- 
agation in the z-direction, both modes turn out to 
have V — ^1/3. For propagation in other directions, 
there are two phonon modes with differing velocities. 



D. Shear modulus for the 2Cube45z crystal 

Let us orient the coordinate axes such 
that {({2} contains the eight wave vec- 
tors (l/-\/3)(±l, ±1, ±1) and {(13} contains 
(l/^/3)(±^/2,0,±l) U (l/%/3)(0, ±v^, ±1). These 
wave vectors are shown in the right panel of Fig. 1. 
The potential energy is given by 
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from which one can read off the nonzero entries of 
the A™""" tensors for / = 2 and / = 3. In the case of 
A2, where the axes are oriented perpendicular to the 
nodal planes of the crystal, the form of A™""" for the 
cubic crystal are easily inferred from the symmetries 
of the cube [6l| . There are in general only three in- 
dependent nonzero entries in A™""", corresponding 
to the terms read from (fTOS]) with the form A^"™™, 
A^""" and A^"™". The form of A J™"" read from 
(|102p is therefore valid to all orders in A2, although 
of course the values of the coefficients, including 
in particular the equality between the A™™"" and 



\ murau 
^2 



obtained from A™""" 
about the z-axis. Note that Ag''^*' 



coefficients, will receive corrections at higher 
order. Finally, note that Ag"""" read from ^02]) is 
by rotating this tensor by 45° 

This is a consequence (after the 45° rotation) of the 
equaHty of A^™"" and A^"™" in the A2 crystal, and 
is therefore not expected to persist at higher order 
in A3. 

As in the previous subsection, we extract the ct/ 
tensors and the matrices of shear moduli vi from the 
potential U of (|102|) , obtaining in this case 





(103) 

As discussed above, there is no symmetry reason for 
= 0, so we expect that this component of the 
shear modulus is nonzero at order A^/(5/i^. Note 
that if we were to rotate our coordinate axes by 45° 
about the z-axis, it would be ^3 that has the zero 
entry while V2 would have all off-diagonal entries 
nonzero. (To confirm this, rotate the A™""" ten- 
sors and re-extract the vi matrices, which are not 
tensors.) 

The A crystal has nodes along x = ((2n + 
l)7r^/3)/(4g), y = ((2?! + l)7rV3)/(4g) and z = 
{{2n + l)7rV3)/(4q). The A3 crystal has nodes 
along x±y ^ ((2n -h l)7rV6)/(4g) and z = {{2n + 
l)7r-\/3)/(4(7). The nodes common to both lie along 
the z = ((2n-|-l)7r-\/3)/(4(7) planes. We therefore ex- 
pect that the crystal will orient itself relative to the 
rotation axis such that rotation vortices lie within 
these planes. Depending on their orientation within 
the planes, they could be pinned where the perpen- 
dicular nodal planes of either the A2 or the A3 crys- 
tals intersect the z = {{2n + 1)11^/3) /{4q) planes. 

We learn from our analysis that the crystals are 
weaker (smaller shear modulus) with respect to 



shear in certain planes. We saw this explicitly for 
fg^, which is zero to order A§ and thus presumably 
weaker although nonzero when higher order terms 
are included. The same will apply to shear in any 
plane obtained from this one by a symmetry trans- 
formation of the crystal, and will apply to the anal- 
ogous planes for the A2 crystal. Note, however, that 
in the 2Cube45z structure the weak planes for the 
A2 and A3 crystals do not coincide. This means 
that if it so happens that motion of a rotational vor- 
tex in a certain direction is only impeded by the 
weaker shear modulus of the A2 crystal, it will in 
fact be obstructed by the stronger shear modulus 
of the A3 crystal, or vice versa. Thus, the relevant 
shear modulus in the analysis of vortex pinning and 
pulsar glitches is the stronger one, which wc find to 
be 16k/9 to order A^. 

As in Section IV. C, we can find further evidence 
for the rigidity of the 2Cube45z crystal by evalu- 
ating the velocity of the transverse and longitudinal 
phonons. Considering the A2-crystal as an example, 
from the dispersion relations ([^^ and P2 = P2 ~ 
P2 = 8k/3 wc find that for propagation in the x- or 
y- or z-dircction the longitudinal phonon mode and 
the two transverse phonon modes all have v = 
For propagation in the x ± z directions, the longitu- 
dinal mode has v = \/2/3 while one transverse mode 

has V = x/T/S and the other transverse mode, cor- 
responding to transverse oscillations for which the 
restoring force would be given by the component of 
the shear modulus which vanishes at order A^ , turns 
out indeed to have v = 0. We see that the velocity 
of both the longitudinal and transverse phonons is 
anisotropic, as expected in a crystal, and see that 
they are comparable in magnitude, confirming that 
the nonzero components of the shear moduli are as 
large as the longitudinal elastic moduli, as expected 
for a very rigid body. 



V. CONCLUSION 

A. The rigidity of crystalline color 
superconducting quark matter 

We have calculated the shear moduli of crys- 
talline color superconducting quark matter with the 
CubcX and 2Cube45z crystal structures. Within 
the Ginzburg-Landau analysis of Ref. [2^, one or 
other of these crystal structures is favored over 
unpaired quark matter and over spatially uniform 
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paired phases like the CFL phase in the wide regime 
of densities given in Eq. ([1]). As we have explained 
in Sections IV. C and IV. D, in both these structures 
the components of the shear moduli that make the 
crystals rigid with respect to vortices pinned within 
them take on the same value to order A-^, given by 



16 

with K defined by ()95|) . Evaluating k yields 



l^CQM 



(104) 



2.47MeV/fm3 ( ^ V f V 



.400 MeV/ 

(105) 

for the shear moduli of crystalline quark matter with 
these two crystal structures. If quark matter is found 
within neutron stars, it is reasonable to estimate 
that its quark number chemical potential will lie in 
the range 



350 MeV < < 500 MeV 



(106) 



The gap parameter A is less well known. Accord- 
ing to the Ginzburg-Landau calculations of Ref. [2^ , 
A/Ao is about 1/4 to 1/2, with Aq the CFL gap 
parameter for Ms = 0. Here, A/Aq for the CubeX 
crystal structure somewhat larger than that for the 
2Cube45z structure and A/Aq is a slowly increas- 
ing function of M^//i, meaning a slowly decreasing 
function of density [20[ . It is reasonable to estimate 
that Ao is between 10 and 100 MeV, but if Ao is 
in the upper half of this range then quark matter at 
accessible densities is likely in the CFL phase, rather 
than in the crystalline phase. So, we suggest that in 
interpreting (jlOSp it is reasonable to estimate that 



5 MeV < A < 25 MeV 



(107) 



keeping in mind that a part of the uncertainty en- 
compassed by this range comes from our lack of 
knowledge of Aq and a part comes from the / fi- 
dependence of A/Aq described in Ref. [13 . The es- 
timates (|107|) and (|106p mean that our result (|105|1 
implies 

0.47 McY/frn^ < i/cqm < 24 MeV/fm^ . (108) 

We shall take this as an estimate of the magnitude 
of i^cQMj although (|105p is a better representation 
of our result for use in future work. 

One qualitative way to appreciate how rigid the 
crystalline phases of quark matter are is to calculate 



the (anisotropic) velocities of long wavelength trans- 
verse and longitudinal phonons, as we have done for 
a few directions of propagation in the CubeX and 
2Cube45z crystal structures in Sections IV C and IV 
D respectively. We find that the transverse modes, 
whose restoring forces are governed by the shear 
moduli, propagate with velocities that are compa- 
rable to the velocity of longitudinal phonons. 

To appreciate more quantitatively how rigid the 
crystalline phases of quark matter prove to be, we 
compare the shear modulus that we have calculated 
to that for the standard neutron star crust, which 
is a conventional crystal of positively charged ions 
immersed in a fluid of electrons (and, at sufficient 
depth, a fluid of neutrons). The shear modulus of 
this solid can be expressed as [6^ 



nj(Ze)^ 



(109) 



where is the number density of ions in the crust, Z 
is the atomic number of the positively charged ions, 
a ~ (3/(47r7ii))^/'^ is the average inter-ion spacing, 
~ 47r/137 and c ^ 0.1 — 0.2 is a dimensionless con- 
stant. Because the crust is electrically neutral, the 
number density of ions is related to rie, the electron 
number density, by rii = n^jZ . And, rie is given in 
terms of the mass and electric chemical potential /le 
of the electrons by 



Tie 



37r2 



(110) 



where fi^ is estimated to be in the range 20 — 80 MeV 
and Z 40 — 50 [fi^l. Using these estimates, we find 



0.092 kcV/fm^ < i^NM < 23 keV/fm^ 



(111) 



Comparing to ()108p . we see that crystalline quark 
matter is more rigid than the conventional neutron 
star crust by at least a factor of 20, and possibly by 
about three orders of magnitude. 

We conclude that crystalline color superconduct- 
ing quark matter is a very good solid indeed, which 
is remarkable since it is at the same time superfiuid. 



B. Toward pulsar glitch phenomenology 

As discussed in the Introduction, the glitches that 
are observed to interrupt the gradual spin-down of 
spinning neutron stars are thought to arise from the 
sudden unpinning of an array of rotational vortices 
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that had been pinned in place, at a fixed area den- 
sity and hence a fixed angular momentum, while 
the other components of the star and in particu- 
lar the observed surface had been gradually slow- 
ing down. When the stressed vortices unpin and 
separate, the superfluid component loses angular 
momentum while the surface spins up. Can these 
phenomena originate within crystalline color super- 
conducting quark matter in the core of a neutron 
star? This phase of matter is a superfluid while 
at the same time having a rigid spatial modulation 
of its superfluid condensate, as we have seen. Un- 
derstanding whether this makes it a plausible locus 
for the origin of pulsar glitches requires addressing 
three questions: Is crystalline color superconduct- 
ing quark matter rigid enough? Do vortices in this 
phase of matter get pinned? And, how rapidly can 
angular momentum be transferred from a crystalline 
quark matter core that has just glitched to the outer 
crust whose surface is observed? 

Our calculation constitutes an affirmative answer 
to the first question. We have shown that both the 
CubcX and 2Cube45z crystal structures have shear 
moduli with magnitude (|105p which are 20 to 1000 
times greater than those of the conventional neutron 
star crust within which glitches have long been as- 
sumed to originate. 

Next, do vortices in fact get pinned? With what 
pinning force? This is a much harder question to 
address quantitatively because doing so requires go- 
ing beyond the long wavelength phonon effective ac- 
tion. The question is what is the difference in the en- 
ergy per unit length of a vortex centered on a nodal 
plane (or at the intersection of two nodal planes) of 
the condensate and one centered half way between 
neighboring nodal planes. Understanding this quan- 
titatively requires constructing a vortex solution in 
the crystalline background, which is a challenging 
task. In the conventional neutron star crust, a vor- 
tex in a neutron superfluid is pinned on "impuri- 
ties" embedded in the superfluid, namely the lat- 
tice of positively charged nuclei. In rotating crys- 
talline color superconducting quark matter, the vor- 
tices are deformations of the phase and magnitude 
of the same condensate whose underlying magnitude 
modulation is the origin of the pinning. Unlike in the 
case of the shear modulus, which describes the re- 
sponse to a stress on length scales long compared to 
those characteristic of the crystal itself, the defor- 
mations introduced by a vortex will occur on length 
scales comparable to the lattice spacing of the un- 



derlying crystal. This means that constructing the 
vortices must be done self-consistently with analyz- 
ing the crystal structure itself — the pinning sites 
are in no sense extraneous impurities. We can pro- 
vide a crude estimate of the pinning force, but we 
defer a quantitative response to this challenge to fu- 
ture work. 

To estimate the pinning force, let us suppose 
(contrafactually) that the core radius of a vortex 
^ ^ 1/A is much smaller than the spacing between 
nodal planes of the crystalline condensate. If such 
a vortex is located where the underlying crystalline 
condensate is maximal, it will have to deform that 
condensate maximally since at the center of the vor- 
tex the condensate must vanish. Clearly, it will be 
energetically advantageous to locate the vortex at 
the intersection of nodal planes where the conden- 
sate already vanishes in the absence of a vortex. This 
argument translates into a pinning energy per unit 
length given at the level of dimensional analysis by 

E 

— ~ f |f^crystallinc| £. (112) 

where |ricrystaiiinc| IS the condensation energy of the 
crystalline phase and where / is some dimensionless 
factor. The corresponding pinning force per unit 
length is given by 

/ l^^crystallincl C /I 1 Q^ 

T " b ' ^^^-^^ 

where the length scale b is half the spacing be- 
tween neighboring nodal planes and hence one quar- 
ter of the lattice spacing. In both the CubeX and 
2Cube45z crystals, b = 7rV3/(4g) = 1.13/S^i. Re- 
calHng that = M^/{8^), we can get a sense 
of the scale of b by seeing that b = 18, 12 fm 
for M'^/fi 100,150 MeV. Reading from plots in 
Ref. [20|, we see that for Aq = 25 MeV this range of 
M'^/fj, corresponds to a robust crystalline phase with 
If^crystaiiinel ~ 2 X 10^ MeV and A - 5 - 10 MeV 
if the crystal has the 2Cube45z structure or A ~ 
10 — 15 MeV if the crystal has the CubeX struc- 
ture. We immediately see that ^ = 1/A and b are 
comparable length scales, which makes this analy- 
sis unreliable at a quantitative level. One way of 
saying this is that the dimensionless factor / must 
then be very much less than one, since the energy 
benefit by moving the vortex by a distance b is of or- 
der |ricrystaiiinc|C^ Only if making this move shifts the 
core from a place where the condensate was maximal 
within the core area to a place where it is close 
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to vanishing. A calculation of / in the case where 
^ ~ 6 as is relevant in our context requires a quanti- 
tative analysis, but it is clear that the energy benefit 
of moving the vortex by a distance b must then be 



J 



^ |f^crystaiiino|?^- Putting the pieces together, we 
can write an estimate of the pinning force per unit 
length as 



0.7 



McV 



(10fm)2 \0.01 



f 



I ^crystalline | 

2 X IQS MeV" 
I 



20 fm 



15 fm 



(114) 



where our choice of / ~ 0.01 as a fiducial value is 
a pure guess and the dependence of the other quan- 
tities in the estimate (|114p on Aq and M^/fi can 
be obtained from the results of Ref. [2^, with the 
fiducial values we have used being reasonable for 
Ao = 25 MeV and M^/i^i = 100 - 150 MeV. 

We can compare our estimate (|114p to the pinning 
force on neutron vortices in a conventional crust j63| , 
in which neutron superfluid vortices are pinned on 

nuclei m, M M, m m, M M m, m, m, m, m, 

[s?! . The pinning energy of a vortex per ion on which 
it is pinned is E^^ - 1-3 MeV BMIil, the ions 
are spaced by a lattice spacing 6nm ~ 25 — 50 fm [i^ , 
and the su per fluid vortices have core radii ^nm ~ 
4 — 20 fm [48| . Hence, the pinning force per unit 
length is [43, 13 



E. 



NM 



1 - 3 MeV 



5nmCnM 



(25 - 50 fm)(4-20fm)' 



(115) 



Although our estimate (|114p is quite uncertain, 
given that we have not constructed the vortex solu- 
tions for rotating crystalline color superconducting 
quark matter, it seems reasonable to estimate that 
the pinning force per unit length on vortices within a 
putative crystalline quark matter neutron star core 
is comparable to that on neutron superfluid vortices 
within a conventional neutron star crust. 

Recent calculations of the profile of vortices in 
BCS-paired superfluid gases of ultracold fermionic 
atoms [g^I may make it easier to estimate the pin- 
ning force on vortices in the crystalline quark matter 
phase. In the cold atom context, it turns out that 
the radius of the vortex core is much smaller than the 
correlation length ^ ~ 1/A which controls the long 
distance form of the vortex profile. Instead, the vor- 
tex core radius is ~ 1/kp, controlled by the Fermi 
momentum rather than by the correlation length. 
If this result were to be obtained in our context, it 



could mean replacing ^ in (|114p by (l//i) ^ 0.5 fm, 
reducing our fiducial estimate (|114p by a factor of 
1600. However, if the vortex cores do turn out to be 
as narrow as this then the assumption with which 
we began our estimate, namely that the core size is 
much less than the lattice spacing, becomes factual 
rather than contrafactual. This would mean that 
there is no longer any reason to expect the dimen- 
sionless factor / to be much smaller than 1, and 
would considerably reduce the uncertainty in the es- 
timate. Replacing / by 1 would increase the esti- 
mate (|114|) by a factor of 100, resulting in a pinning 
force which is only slightly smaller that on the super- 
fluid neutron vortices within a conventional neutron 
star crust. We leave the determination of the pro- 
file of vortices in crystalline quark matter to future 
work, but it will clearly be very interesting to see 
whether they have narrow cores as in Ref. [6J|, and 
if so whether their pinning turns out to be controlled 
by their core radii or by the correlation length. 

The third question which must be addressed is 
how, and how quickly, angular momentum can be 
transferred from a crystalline quark matter core to 
the outer crust. Some glitches arc known to oc- 
cur on timcscalcs of minutes which means that if 
a glitch occurs within the core angular momentum 
must be transferred to the observed crust at least 
this fast. The core and crust are linked via being 
bathed in the same electron fluid and via magnetic 
fields. In the conventional glitch scenario, when the 
neutron superfluid in the crust suddenly slows down 
as its vortices come unpinned and the nonsuperfluid 
component of the crust, which includes the ions and 
the electrons, speeds up, the electron fluid couples 
the crust to the core well enough that the core also 
speeds up within seconds [i^, [gBl • We therefore ex- 
pect that if a glitch occurs within a crystalline quark 
matter core, with this superfluid component slow- 
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ing down, and if moving vortices can impart angular 
momentum to the electrons, then the electron fluid 
will ensure that the entire rest of the star including 
the outer crust speeds up. In the conventional sce- 
nario, the mechanism by which moving vortices exert 
a torque on the ions, and hence the electrons, in the 
crust has been described in Refs. [5l|, [H, [s^l • In our 
case, we have not demonstrated how moving vor- 
tices in the crystalline phase can torque up the fluid 
of gapless charged quark quasiparticles, and hence 
the electrons. 

So, tallying the status of the three questions that 
must be addressed: The first is settled, answered 
in the affirmative by our calculation of the shear 
modulus of the crystalline phase. The second re- 
mains to be addressed quantitatively but seems to 
be answered in the affirmative by our crude esti- 
mate of the pinning force on vortices in the crys- 
talline phase. The third remains open, not yet ad- 
dressed in a satisfactory fashion but nevertheless 
with no reason to doubt that its answer is also affir- 
mative. Addressing the third question and address- 
ing the second question quantitatively both require 
constructing the rotational vortex solutions for ro- 
tating crystalline quark matter. This is therefore the 
crucial remaining step in completing the connection 
between the microphysics of crystalline color super- 
conducting quark matter and the phenomenology of 
pulsar glitches, and hence determining whether the 
characteristics of observed glitches rule out, or are 
consistent with, the presence of crystalline color su- 
perconducting quark matter within neutron stars. 

It is also worth asking whether a "core-quake" 
scenario could be a viable model of glitches, qual- 
itatively distinct from that based upon vortex pin- 
ning [40 . 66„ 6^ . As a spinning neutron star slows 
down, it becomes less oblate. This will require 
macroscopic adjustments to the shape of a putative 
crystalline quark matter core. Given the enormous 
shear moduli of this rigid phase of matter, enormous 
amounts of elastic energy would be stored as the core 
is deformed and stressed, energy which would be re- 
leased in core-quakes during which the crystalline 
core "breaks" and rearranges its structure so as to 
reduce its moment of inertia, consequently increas- 
ing its angular velocity. The original "crust-quake" 
model for pulsar glitches failed because it failed 
to describe the magnitude and frequency of glitches 
in the Vela pulsar [H, [13, [H, [53|. Now that we 
know that crystalline quark matter has shear mod- 
uli which are 20 to 1000 times larger than those of 



the crust, core-quakes are worth re-investigating. 

Finally, the advent of gravity wave detectors opens 
a new possibility for unique signatures of the pres- 
ence of rigid matter within neutron stars, indepen- 
dent of transient phenomena like glitches [6^ . LIGO 
has already set limits on the steady-state gravity 
wave emission from 78 nearby pulsars (70| , and these 
limits will become more stringent. Owen's work (69j 
shows that if an entire neutron star were solid, with a 
shear modulus as large as that we find for crystalline 
color superconducting quark matter, it could in prin- 
ciple support a quadrupole moment large enough 
such that the resulting gravity waves would already 
have been detected. This suggests that as the ob- 
servational upper limits improve, the size of a pu- 
tative rigid crystalline color superconducting quark 
matter core could be constrained. However, Owen's 
estimates for a star that is rigid in its entirety can- 
not be applied straightforwardly to the case where 
there is a rigid core surrounded by a fluid "mantle" . 
Oblateness about the rotational z-axis is not enough 
to generate gravity waves; the quadrupole moment 
must be mis-aligned, such that the moment of inertia 
tensor satisfies I^x 7^ lyy It is hard to imagine how 
this could come about for a rigid quark matter core, 
whose shape will equilibrate to follow an equipoten- 
tial surface via converting core material into mantle 
material or vice versa as needed at different loca- 
tions along the core/mantle interface. Nevertheless, 
this line of enquiry warrants careful investigation. 
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APPENDIX A: PHONON MASS IS ZERO TO 
ALL ORDERS IN A 

In this Appendix, wc show explicitly that there 
cannot be any term in the phonon effective action 
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at any order in A and at any order in u/(a;) which 
is nonzero if Uj (x) is constant over space and time, 
other than a trivial constant. Equivalently, we show 
that if we expand the Lagrangian density as a power 
series in u/(a;), every term has at least one derivative 
acting on each ui(x), meaning that the mass of the 
phonons is zero. This is guaranteed by Goldstone's 
theorem, but the explicit demonstration of this "ob- 
vious" result is nontrivial and for this reason we give 
it here. 

Before proceeding, we recall Eqs. (^5]) and (^5]) . 
which imply that 

AJix) - ^ cxp(2iq? • (r - ui{x)) . (Al) 
q? 

Since we want to prove the result to arbitrary pow- 
ers in the u fields, we do not make an expansion in 
small u. In terms of Feynman diagrams, this means 
that in evaluating the effective action in Eq. (|43)l 
we resum the vertices with increasing powers of u/. 
Therefore any vertex for A/ comes with a factor 
siapsiij exp(— 2iqj • u/(a;)) and a momentum inser- 
tion 2qf. On the other hand, the vertex for Aj 
comes with a vertex factor e/a/jE/y exp(2iqj • u/(a;)) 



and a momentum insertion — 2qj. 

Integrating out the x fields in Eq. (|^ is equiva- 
lent to calculating all possible one-fermion-loop dia- 
grams with arbitrarily many external phonon inser- 
tions. (Higher loop diagrams with internal phonon 
propagators are suppressed by powers of Air/ A, 
with Air the typical energy of the phonon fields, and 
are therefore completely negligible.) A generic one- 
loop diagram will have nj vertices proportional to 
A/ exp(— 2iqj'' • u/(x)) at which momenta 2qj'' are 
inserted into the loop, and nj vertices proportional 
to A|exp(2iqJ^ • u/(x)) at which momenta — 2qJ^ 
are inserted into the loop. Here, k = 1, 2, 3...nj and 
r = 1,2, 3.. n/. The number of appearances of A/ 
has to be the same as that of AJ because the diagram 
can only depend on powers of |A/p. Different one- 
loop diagrams correspond to different choices of ni, 
71,2 and ^3 and different choices of the q"" 's and the 
qj^ 's from among the sets {q/}. The color and flavor 
indices for the propagators in the diagram linking 
the vertices into a loop are chosen to be consistent 
with the color and flavor epsilon symbols associated 
to each vertex. The contribution from a generic one- 
loop diagram will be 



Xoc |AiP"i|A2p"^|A3p"3 X Tr 



1 ^2q°i.(r-ui(a:)) ^ ^-2q:i ■(r-U2(x)) ^ ^2q°" -(r-ui (a:)) ^ 



(A2) 

We are interested in evaluating only the contribution from such a diagram in which no derivatives act on 
any u/ fields. This contribution is given by 



I(x |Ai|^"i|A2r"^|A3|^"-^ X Tr 



exp{z5^2u,(x).(^q?-^q5)} 



a—ai b—bi 



1 



1 



= -2q2'-r. 



I 



.2qr^r 



i9 + /i.,3 



(A3) 



Momentum conservation implies that the net mo- 
mentum added to the loop is zero, i.e. 

/ \ a=ai b=bi I 

(As we argued after Eq. (|54p . the momenta con- 
tributed by (j) is much smaller than | q^ | . ) 



We now recall from Section II that the magnitude 
of the qj are different for different /. (jqjl = ??(5m2 
is close in value to jqgl = rjSy,^, but they are not 
exactly equal because 5^2 and <5/^3 differ by terms 
of order M^/n^.) This means that the momentum 
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conservation condition (|A4[) can only be satisfied if 
E q? - E = , (A5) 

a—ai h—bi 

separately for each I. This implies that in Eq. (|A3p . 
the coeflacients of each of the U/ cancel, and thus 
implies that Eq. (jA3[) . the contribution of a generic 
one-loop diagram in which no derivatives act on any 
u/'s, is independent of u/, making it a trivial con- 



stant in the phonon effective action. The phonons 
arc therefore massless to all orders in A and u/. 

As a special case, we can now demonstrate 
Eq. ([?T|) explicitly. That is, we can show explic- 
itly that -4^/ d'*'xtTcF{{'!^cF)^ ^cf) is indepen- 
dent of U/ as must be the case since it includes no 
derivatives acting on u/, and so would constitute a 
mass term for the u/ if it were to depend on the u/. 
Indeed, 



^ I d^xti:cF{iA"cF)^A-cF) 



°-(r-u/(a;))g-2qj-(r-u/(a:)) 



(A6) 



r 



as given in ([5 



APPENDIX B: SINGLE PLANE WAVE 

In this Appendix, we investigate phonons in the 
presence of a condensate for which only one of the 
three A/ is nonzero, meaning that only quarks with 
two different colors and flavors form Cooper pairs. 
We take A3 ^ and Ai A2 = 0. This imphes 
that there is pairing only between ur and dg quarks, 
and between ug and dr quarks, where r, g and b refer 
to the colors of the quarks. Furthermore, we assume 
that the set {q^} contains only one vector, q, so 
that A3(x) varies in space as a single plane wave, 
A3(r) = A3exp(2iq • r). In this simple case it is 
possible to derive the phonon effective action and 
the shear moduli without employing the Ginzburg- 
Landau expansion, working to order c/)^. 

The oscillation of the condensate is described by a 
single phonon field, U3(a;). In analogy with Eq. ([26|) . 
we define 

A^{x) = A3e2*i-('"""-^(^» 

« A^e'^'i '^ (^1 - icj^ix) - \m^^ (Bl) 

where (j){x) = 2q • U3(a;). 

This condensate breaks translational symmetry in 
the q direction, but leaves an 0(2) symmetry cor- 



responding to rotations about the q-axis unbroken. 
Taking q along the z axis, the potential energy U of 
Eqs. ((82l) and ((84l) must be symmetric under rota- 
tions about the z-axis, taking the form 



U =- 



(B2) 



where the strain tensors s™" are defined in (|83p . 
y^zzzz yxzxz _ y%jzyz ^.^^ independent 

clastic moduli that we will now evaluate. 

In writing the Lagrangian in terms of the Nambu- 
Gorkov fields, we ignore the quarks that do not par- 
ticipate in pairing. We also note that the inverse 
propagator can be written as a block diagonal ma- 
trix made up of four blocks that correspond to the 
ur particles and dg holes, dg particles and ur holes, 
ug particles and dr holes and dr particles and ug 
holes. Since only quarks that belong to the same 
block interact, the inverse Nambu-Gorkov propaga- 
tor can be written as 



^ ^ ~ ^ur-dg ' 



with 



0-1 

dg—ur 



s: 



(B3) 
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ur — dg 



ug—dr 



' '-^dg-ur 



' dr—uq 



*^ + ^, (C^7')A3"(^) 

(C7^)Ar(x) 



(B4) 



The phonon effective action is obtained by integrat- 
ing out the Nambu-Gorkov fields, yielding the result 



z5[u]=log(Z[u]) = -zi(FT)|A 



2Tr, 



S. 



-1 

ur — dg 



(B5) 



where we have used the property that the contribu- 



tions from the four blocks, Sj_^g, S^^_^^, Sj_^g 

and S^g_^^, are equal. Since the trace of an opera- 
tor does not change upon making a unitary transfor- 
mation, corresponding to a change of basis, we can 
simplify S^^ by choosing a basis which gets rid of 
the exp(2iq • r) and the Cj^ appearing with A3. In 
the new basis, 



m — dg 



e-'i"" 

{C-/^)e" 



A3(l-10^ 



e*'!"" 

(-C75)e-"i-'' 



(B6) 



In this simple case, the inverse propagator in the 
absence of the phonons can be inverted. (This is 
why we do not need to resort to the Ginzburg Lan- 
dau approximation.) Hence we separate S~^_, into 



-dg 



c-l 

'-'a 



S with 



and 

E = 



A5(+*0-i</>2) 



A3 



A3H</,-i02) 



(B7) 



(B8) 



Since we have rotated out the phase exp(2iq-r) from 
A3, the off-diagonal components of do not de- 
pend on position anymore. Upon inverting 5^^, one 
gets the full propagator for the fermions in the ab- 
sence of phonons, to all orders in A. This propa- 
gator is diagonal in momentum space and therefore 
can be written in a simple way employing the weak- 
coupling approximation and the HDET formalism 
(implemented as in Eqs. ([SS]) and (|5D|) or more for- 
mally as in Ref. [lOl), yielding 



S'a(p) 



DAip) 



(F • p - V • q - (5/i3) -A3 

-A^ (F-p-v.q-<5Ai3) 



(B9) 



where 

DAip) 



(y • p - V • q - Sfi3){V • p - V • q - Sfia) 
-lAal^ (BIO) 
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a) 



p+k 



b) 



FIG. 5: Diagrams contributing to the phonon self-energy. The dot-dashed hnes correspond to the propagator of 
the phonon fields. Full lines correspond to the full propagator of quark quasi-particles. The phonon-quark-quark 
vertices Fi and the phonon-phonon-quark-quark vertex F2 are also shown. Diagram b) contributes only to the zero 
momentum polarization tensor. Diagram a) is nonvanishing for any value of the external momentum k of the phonon. 



We now consider the interaction with the phonon 
field. To second order in 0, there are two interaction 
vertices, i(j){x)ri and ~^(j){x)'^T2, where 

(BU) 

The phonon effective action to order (j)'^ is given by 



J 



the diagrams in Fig. [5l namely 

5*' = Tr(5A ^Ti^A 0ri) - Tr(5A Tz) ,(B12) 

where the trace is over Nambu-Gorkov indices and 
space-time. Evaluating the action along the same 
lines as in Section UlI El wc find 



dv {V-k){V-k) 
1^ DAip + k)DAip) 



(B13) 



where the fc- independent second diagram in Fig.[5]has cancelled the fc-independent contribution from the first 
diagram. We then integrate over pq and v analytically and do the s-integral numerically. Upon returning to 
position space, the action takes the form (|90p which, in this simple case, reduces to 



4g' 



(B14) 



where 9|| = q(q • d) and d± = d — We have 
written the factors of = Arf'SiJ?' in the denomi- 
nators in (|B14p in order to facilitate comparison to 
(|B2p . The potential energy U in (|B14p can easily be 
written in the form given in Eq. (|B2p by substituting 
(/) = 2q • U3 and keeping only the term symmetric in 
diXXj. Taking q in the z direction, one can easily see 
that A^^^^ = A|| and A"=^"=^ = A^^w^ = A_l. 

The coefficients p, Ax and A|| are each propor- 
tional to //^ and have nontrivial A3- and (5/i3- 



dependence. We have obtained them numerically for 
arbitrary values of A3, and by plotting them versus 
A§ at small A3 we have checked that p and A|| are 
proportional to A3 and A_l is consistent with be- 
ing proportional to A|. Within the accuracy of our 
numerical analysis, we find p = A|| — k to order 
A3, where k is given by ([M]) . This is in agreement 
with what we obtain from the phonon effective ac- 
tion evaluated to order A| in the Ginzburg- 
Landau approximation, upon specializing (j93[) to 
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the single plane wave case. It therefore provides a 
useful check on our calculations. Whether via the 
Ginzburg-Landau analysis of the main Sections of 
this paper or via the numerical all-order-in-As eval- 
uation of this Appendix, we learn that A^^^^ = k and 



= to order A3, from which we conclude that 
the shear modulus of this singic-plane-wave "crys- 
tal" is zero, to this order, and is presumably propor- 
tional to ^^A|/(S/i^. 
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